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The present contribution is aimed at removing most of the obstacles in understanding the quan- 
tum mechanics of massive particles. We advance the opinion that the probabilistic character of 
quantum mechanics does not originate from uncertainties caused by the process of measurement 
or observation, but rather reflects the presence of objectively existing vacuum fluctuations whose 
action on massive particles is calibrated by Planck's constant and effects an additional irregular 
motion. As in the theory of diffusion the behavior of a single particle will be described by an 
ensemble of identically prepared but statistically independent one-particle systems. Energy con- 
servation despite the occurrence of a Brownian-type additional motion is achieved by subdividing 
the ensemble into two equally large sub-ensembles for each of which one obtains an equation of 
motion that has the form of a Navier-Stokes- or "anti"-Navier-Stokes-type equation, respectively. 
By averaging over the total ensemble one obtains a new equation of motion which can be converted 
into the time-dependent Schrodinger equation. We clarify the problem of the uniqueness of the wave 
function and the quantization of orbital momentum. The concept allows the inclusion of electro- 
magnetic fields and can be extended to interacting N-particle systems. We analyze the problem of 
how an experimental setup can consistently be decomposed into the quantum system under study 
and the residual quantum system "apparatus". The irregular extra motion of the particle under 
study allows a decomposition of the associated ensemble into two subensembles the members of 
which perform, respectively, a right-handed or left-handed irregular circular motion about a given 
axis which becomes physically relevant in the presence of a magnetic field. We demonstrate that 
this orientation-decomposed "Zitterbewegung" behaves - in accordance with Schrodinger's original 
idea - as a spin-type angular momentum which appears in addition to a possible orbital angular 
moment of the particle. We derive the non-relativistic time-dependent Pauli equation and propose 
a theory of the Stern-Gerlach experiment. The Dirac equation proves to be derivable by drawing on 
similar arguments used in obtaining the Pauli equation. We, further, attempt to put Bell's theorem 
and the Kochen-Specker theorem into perspective. 
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"... to skeptics, heretics and naive realists everywhere. 
Keep doubting; let others keep the faith." 

David Wick in: The Infamous Boundary 1 

"I have never been able to discover any well-founded rea- 
sons as to why there exists so high a degree of confidence 
in the current form of quantum theory." 

David Bohm in: Wholeness and the Implicate Order 2 



I. INTRODUCTION 

There exists a rich literature on attempts that have 
been made to derive non-relativistic quantum mechan- 
ics from a concept of dissipationless stochastic point me- 
chanics. A precursor of the idea of correlating the prob- 
abilistic character of quantum mechanics with the action 
of a stochastic background field may be seen in the paper 
by Bohm and Vigier— . The present contribution draws 
on later work on this subject but avoids certain impli- 
cations that have often been criticized during the past 
20 years. (S. e.g. W. Weizcl 4 , E. Nelson^, Guerra and 
Morato, 7 , M. BaublitaS L. de la Pena and A. H. Cetto^, 
Petroni and Morato™ T. C. Wallstromii and numerous 
references therein. For a rather complete review see R.F. 
Streater—. An earlier review covering work up to 1986 is 
given in a book by Namsrai^.) A more recent contribu- 
tion is due to Fritsche and Haugk^. 
The derivation of the time-dependent Schrodinger equa- 
tion constitutes the focus of the following considera- 
tions. A very interesting alternative to our approach 
that also relates to vacuum fluctuations, but draws on 
non-equilibrium thermodynamics has recently been put 
forward by Grossingi^ 1 ^. Based on the concept of our 
derivation one is led to conclude that every conceivable 
situation of a physical system is exhaustively described 
by the respective solution to the Schrodinger equation, 



and that there can be no independent measurement prob- 
lem. As for this point we side with J. Belli 7 - who argues 
that the attempt to base the interpretation of quantum 
mechanics on some notion of "measurement" has raised 
more problems than it has solved. As we shall outline in 
Section IXXIII "measurements" relate outcomes, e. g. de- 
tector readings, to characteristic properties of a quantum 
system by using solutions to the Schrodinger (or Pauli) 
equation as primordeal information. Without these equa- 
tions and their solutions "measurements" , i. e. in general, 
detector or "pointer" readings, constitute a set of worth- 
less data. 



II. ORIGIN OF QUANTUM MECHANICAL 
RANDOMNESS 

We interpret the fact that microscopic particles move 
and behave differently from macroscopic objects as re- 
flecting the active role of the vacuum providing a back- 
ground for energy fluctuations. The latter will henceforth 
be referred to as vacuum fluctuations. The consequences 
of their existence have already been discussed quite 
some time ago, s. e.g. Bess^, Puthoff 19 , 20 , Boyer 21 , 
Calogero— Carati and Calgani 2 ^. Present day quantum 
mechanics is strongly shaped by historical contingencies 
in its development, and it has become almost impossible 
to tell fiction from facts. Statements on "the measure- 
ment of positions at different times" and "there is no 
momentum of a particle in advance of its measurement" 
are typical of this school of thought (s. e.g. Streater— ), 
yet they are definitely void of meaning. What kind of 
experimental setup should allow a perfectly accurate po- 
sition measurement at a perfectly accurate time point? 
And how does the setup look like that allows the mea- 
surement of a particle momentum in the spirit of ortho- 
dox quantum mechanics; i. e. with zero variance? It 
is totally impossible to perform non-fictional measure- 
ments on quantities that would conform to their quan- 
tum mechanical definition, e. g. measuring commuting 
observables like energy and angular momentum at the 
same time. The "observables" around which a substan- 
tial portion of quantum mechanical literature revolves are 
in reality non-observables. Further, there is simply no ev- 
idence of a causal interrelation between the probabilistic 
character of quantum mechanics and indeterminacies in- 
troduced by "the observer" . 

By contrast, there is every reason to believe that vacuum 
fluctuations are real and constitute an objective prop- 
erty of nature. Zero point motion of particles constitutes 
the most obvious evidence of their existence. It is this 
zero point motion which, for example, keeps liquid 4 He 
"molten" down to the very lowest temperatures and ex- 
plains this extraordinary material property. 
One could view vacuum fluctuations as caused by an ex- 
change of energy between the mechanical system in ques- 
tion and the embedding vacuum that serves as an energy 
reservoir in terms of virtual particles: if that reservoir 
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reduces its content of virtual particles, the energy of the 
system under study increases so that the energy of the 
entire system comprising this "vacuum reservoir" is con- 
served. Considerations of Calogerc 2 — point in a similar 
direction. In that sense quantum mechanical systems 
may be viewed as open systems like classical point mass 
systems in contact with a heat bath. This analogy will 
become particularly visible in our treatment. We shall 
use the terms "point mass" and "point charge" with the 
reservation that the actual size of the particles in question 
might well be finite of the order 10 -13 cm, but very small 
compared to atomic diameters of the order 10~ 8 cm. Oc- 
casionally "point mass" will stand for the centroid of an 
atom or some composite system. 

In the following we shall focus on the description of the 
subsystem "point mass in real-space" which is open to- 
ward the active vacuum and whose energy is therefore 
conserved only on average. 

An implication of this concept is that charged point 
masses, despite their irregular motion, do not emit or 
absorb radiation on the average. In stationary situations 
a charged point mass will exchange photons with the vac- 
uum in a way that does not change its average energy and 
momentum. 

Radiation only occurs when the probability density of the 
point mass, being at its various positions in space, or the 
associated current density becomes time-dependent. 

This is analogous to a system kept by non-heat conduct- 
ing fibers in a vacuum chamber whose walls serve as a 
heat bath. In a stationary state situation the system ex- 
changes constantly photons with the heat bath without 
changing its average energy. However, if its temperature 
is, for example, higher than that of the wall, the system 
starts radiating, that is, there is now a net flow of pho- 
tons leaving the system. 

If one disregards the details of the energy transfer be- 
tween the two systems, vacuum fluctuations appear as an 
irregular temporary departure of the particle in question 
from its energy conserving trajectory in that it changes 
its energy by an average amount AE for an average time 
interval At so that AE At — f ft where h = 2 irh denotes 
Planck's constant and the factor / is about unity. It is 
this departure from classical energy conservation which 
explains, as already alluded to, why a harmonic oscillator 
in its state of lowest energy is irregularly driven out of the 
position where it would be classically at rest. Further- 
more, it explains the stability of a hydrogen atom in its 
ground-state (which applies quite generally to all atoms 
and their compounds), the zero-point motion of atoms 
in molecules and solids and the "tunneling" of particles 
through a potential wall which actually amounts to over- 
coming that wall. 

Zero-point motion is commonly associated with the un- 
certainty relation which, however, merely shifts the prob- 
lem of understanding a non-classical phenomenon to un- 
derstanding the origin of a non-classical relation. More- 



over, it amounts to keeping a blind eye on the fact 
that one is dealing here with a ground-state phenomenon 
which is certainly not observer-induced. Only if the con- 
trary would apply, one would be justified in referring to 
the uncertainty relation. 

A Boltzmann distribution over the energy levels of some 
system is completely independent of the details and the 
kind of the energy exchange between the heat bath and 
the system. The distribution contains only one univer- 
sal parameter, viz. Boltzmann's constant. Similarly, a 
system's stationary zero-temperature states that emerge 
from exchanging energy with the vacuum do not depend 
on the details of this exchange and on the kind of particles 
involved, but only depend on another universal constant, 
viz. Planck's constant. 

The envisaged derivation implies that particle trajecto- 
ries persist under the influence of the stochastic vacuum 
forces. Their existence becomes particularly obvious with 
tracks of a-particles in a track chamber, but also with the 
trajectories of electrons in a field electron microscope. 
Their property of forming straight lines from the field- 
emission tip (assumed semi-spherical) to the monitoring 
screen is actually presupposed in calculating the magni- 
fication of the microscope. Conversely, purely quantum 
mechanical behavior occurs at lowest energies when the 
trajectories do no longer possess a classical reference in 
the limit % — > 0. Trajectories still persist in that case, 
but the respective particle now performs a purely irreg- 
ular motion. 

The existence of particle trajectories is denied by the 
Copenhagen school of thought because "things that can- 
not be observed do not exist" . Supporter of this view 
have to live with the conflict that a complex- valued wave- 
function or its associated state vector, which constitutes 
the center of quantum mechanics, cannot be observed as 
well. By contrast, we believe that the validity of assump- 
tions can only be scrutinized by checking the consistency 
of the resulting theory against experimentally accessible 
quantities and laws. We are here in complete accord with 
Ballentine who states in his seminal article^: 

"...quantum theory is not inconsistent with the suppo- 
sition that a particle has at any instant both a definite 
position and a definite momentum, although there is a 
widespread folklore to the contrary. " 

In Section IIIII we briefly discuss the construction of en- 
semble averages of quantities that appear in the Navier- 
Stokes equation given in Section[lV] We regard this equa- 
tion as a mathematical object that derives entirely from 
classical concepts. Details of its derivation, which goes 
essentially back to Gebelein^, will be relegated to the 
Appendix, Section fXXXIIl We discuss the construction 
of a "Brownian" and an "anti-Brownian" sub-ensemble. 
The motional behavior of the latter is governed by an 
"anti-Navier-Stokes" equation. We explain why a system 
of statistically independent particles moves according to 
the arithmetic mean of these two equations when their 
motion is governed by classical mechanics plus "conserva- 
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tive" stochastic forces. On forming this arithmetic mean 
we arrive at an equation that can be converted into the 
Schrodinger equation. We demonstrate that Wallstrom's 
objection! 1 against the legitimacy of this conversion and 
his arguments in favor of the standard approach to the 
quantization of orbital momentum are based on a mis- 
understanding and ignore fundamental considerations of 
Pauli2& and Born and Jordan^! in the early days of "con- 
ventional" quantum mechanics. In Section IXIVI we show 
how the derivation of the time-dependent Schrodinger 
equation can be extended by including electromagnetic 
fields. The derivation can be extended further to inter- 
acting many-particle systems. 



III. DEFINING ENSEMBLES AND AVERAGES 

As in the theory of diffusion we start with considering 
a point-like particle that is driven by an external con- 
servative force F(r) and moves in an environment where 
it is exposed to additional stochastic forces. To gain ac- 
cess to quantities that are commonly discussed within 
this framework we construct a sufficiently large set of N 
identical systems (an ensemble of systems) under the sup- 
position that there is no correlation between the stochas- 
tic forces of different systems. As a fundamental con- 
sequence, one is led then, as in the theory of diffusion, 
to a form of quantum mechanics that merely describes 
ensemble behavior. But this, again, is in accord with 
Ballentine's viewer "..in general, quantum theory pre- 
dicts nothing which is relevant to a single measurement 
( excluding strict conservation laws like those of charge, 
energy or momentum)." 

The relative freqency with which the particle appears at 
the time t in an elementary volume A 3 r around the point 
r is given by 



n(r, t) 
N 



p(r,t) A 3 r 



(1) 



where n(r, t) is the number of particles in A 3 r, and 
p(r, t) denotes the probability density. We, furthermore, 
introduce N r for the number of elementary volumes into 
which the total volume V is thought to be subdivided. 
Since the sum over all elementary cells yields N particles 
we have 



En(r.t) 
N 



1 that is / p(r,t)d 3 r = 1 
Jv 



(2) 



We refrain here from discussing the proper limiting case 
N — > oo and relating relative frequencies to probabilities, 
as this matter has extensively been analyzed elsewhere (s. 
e. g. Streater— ). We assume that there will always be a 
smooth function p{r,t) for any finite N that provides 
a least mean square fit to the actually histogram-type 
function in real-space where j numbers the cubes 

into which the normalization volume V is thought to be 



subdivided, and Tj denotes the centroid of the particle 
positions in the respective cube. 



The relative frequency 



%(r,t) 
N 



which we shall below ex- 



press as the mod squared of some wave function ij){r,t), 
refers - when multiplied by A 3 r - to the subset of iden- 
tically prepared systems where the particle appears at 
r and nowhere else simultaneously, otherwise the term 
"particle" would be meaningless. We think that the com- 
monly used phraseology "probability of finding the par- 
ticle at r" is inappropriate because it suggests that one 
would have placed a detector at r monitoring the occur- 
rence of that particle. However, a detector would - apart 
from causing various uncontrollable perturbations - ter- 
minate the motion of the particle on impact, and hence 
there would be a shadow area behind the detector where 
n ( r,t > 0, different from the original unperturbed situ- 



ation. Wherever in the following the quantity 



i(r,t) 
N 



A 3 r 



or p(r,t) A 3 r will appear it is clearly to be understood 
as the probability of the particle being in A 3 r around r. 

We temporarily number the particles in A 3 r at time t 
by an index i, (i = 1, 2 . . . n(r, t)). The particles move, 
in general, at different velocities Vi(t). We define the 
ensemble average of the latter as 



v(r,t) 



1 



i(r,t) 



n(r,t) ^ 



(3) 



As is familiar from the theory of diffusion, the individ- 
ual velocities Vi(t) in A 3 r will in general be quite dif- 
ferent from v(r,t) which we shall come back to later in 
Section [Xj By contrast, in Bohm's version of quantum 
mechanics^ the true particle trajectories are, for no ob- 
vious reason, identified with the streamlines of the ve- 
locity field v(r,t). This is one of the points where our 
approach differs fundamentally from Bohm's and reflects 
a concomitant feature of our definition of v(r,t): 
In performing the average according to Eq.© one sums 
over velocities Vi(t) of different trajectories that run 
through sometimes very different regions of the avail- 
able space of the one-particle system. Hence, they are 
influenced by the classical field F(r) in those regions. 
This carries over to the ensemble average v(r,t). That 
means: if one places a diaphragm somewhere so that a 
continuous subset of trajectories is blocked out, v(r,t) 
changes. That kind of non-local sensitivity explains why 
the streamlines of the field v(r, t) are affected by portions 
of the space which may be far away. The unfamiliar fea- 
ture of non-locality will be illustrated by a particularly 
surprising example in Section [IXI 

As a general property of the stochastic forces that act on 
the respective particle in each system, we require them 
to ensure ergodicity in the following sense: 
If the system is not explicitly time-dependent, that is, 
when it is in a bound stationary state and if one would 
follow the particle on its trajectory within the range it is 
bound to, one would see it successively occur in all the 
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cubes A 3 r over which - in the ensemble average - all par- 
ticles of the ensemble are distributed at a certain instant 
t. Thus, instead of forming the ensemble average accord- 
ing to Eq.([T|) it can for a single particle just as well be 
defined as 



lim 



At(r) 
T 



p{r) A 3 



(4) 



where At(r) denotes the overall time which the particle 
has spent occurring repeatedly in A 3 r around r within 
the total time span T. 

The velocity v(r) can be defined analogously 



v(r) 



1 



h(r) 



n(r) 

E 

i=l 



(5) 



where n(r) is the number of times the particle has 
occurred in A 3 r around r, and U denotes some point 
within the time span the particle has spent there the i th 
time. In realistic cases in which the system under study 
undergoes transitions between quasi-stationary states, 
one has to allow T to be finite, and quasi-stationarity 
can only be ensured if the changes are sufficiently slow 
on a time scale of unit length T. Practical experience 
shows, that this applies to the majority of cases. How- 
ever, in Section IXXII we shall give an example where 
T must be expected to be far too long to justify a 
classification of the states in a photo emission transition 
as quasi-stationary. 

Yet, the bulk of this article will deal with ensemble 
averages. 



IV. NAVIER-STOKES EQUATIONS 

If a particle of mass mo moves in an environment of 
kinematic viscosity v the resulting ensemble average of its 
velocity v(r, t) is just the sum of the so-called "convective 
velocity" v c (r, t) and a "diffusive velocity" u(r, t) driven 
by the stochastic forces of the embedding medium: 



v(r, t) = v c (r, t) + u(r, t) 



(6) 



Employing the Smoluchowski equation (s. Section 
IXXXIII) for the probability density p(r,t), similarly for 
the probability current density j c (r,t) = p(r,t)v c (r,t) 
and invoking Einstein's laws 3 - for the mean square dis- 
placement we obtain a Navier-Stokes-type equation of 
the form 



d_ 

di 



(v — u) + [(v + u) ■ V(u — u)] — v A(y — u) 



too 



F(r). (7) 



with F (r) = — VV r (r) denoting the external conservative 
force acting on the particle. The "osmotic" or "diffusive" 



velocity u(r,t) is defined by 



u(r, t) = —v 



p{r,t) 



(8) 



or equivalently in terms of the diffusion current density 
3d 



where 



j D (r,t) = -vVp{r,t) 'Tick's law" (9) 



j D (r,t)=p(r,t)u(r,t). (10) 



In the special case when v c = the equation of continuity 
reduces to 



0_ 

di 



p + V • pu = 0, 



(11) 



which on insertion of u(r, t) from Eq.© attains the form 
of the diffusion equation 



d_ 

di 



p = v Ap . 



(12) 



On the other hand, when v = one has u(r, t) = 0, and 
hence all particles move now along smooth trajectories 
r(t) so that the various velocities Vi(t) under the sum in 
Eq.Q become equal: Vi(t) = v(r(t)). Thus 

d 

— v(r(t)) = (k = 1,2,3) ->D'V»eO, 

OXk 

and consequently Eq.([7]) reduces to Newton's second law. 

The set of equations ||7J) to ([9]) will be derived in Section 
IXXXIII 

Eq.© may be rewritten 

u(r, t) = - V YP^11 = ^ Hp{r . t )/ Po ] (13) 

where po denotes a constant density that has merely been 
inserted for dimensional reasons. As u(r, t) can be ex- 
pressed as a gradient of a function, we have 



V x u(r,t) = , 



and hence 



(u ■ V)« = V 



u 



(14) 



(15) 



If we, further, make use of the identity 

V x (V x a) = V(V • a) - Aa (16) 

and observe Eq.(fT4|) we obtain Au(r, t) — V(V • u(r, t)). 
Thus, Eq.((7]) in conjunction with Eq.® may be cast as 



TO -v(r, t) = F(r) - VV stoch {r, t) + Q(r, t) , (17) 
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where V r s t oc / l (r, t) and ft(r,t) are abbreviations which 
stand for 



Vstoch — v 



and 



1 fVp\ 2 X7 2 p 

2 ^7 



(it • V) w + v A w . 



(18) 



In deriving (fT5|) we have observed that A Vit 
(^) 2 . Furthermore, we have introduced ^ as the "con- 
vective (or hydrodynamic) acceleration" which in the 
present context merely represents an abbreviation 



dv(r,t) dv 



dt 



dt 



+ v ■ Vt> . 



(19) 



The "stochastic potential" V s toch(r, t) depends on v 2 
whereas Q(r,t) is proportional to v. 
The latter constant is associated with the occurrence of 
the stochastic forces which - in the absence of an exter- 
nal force F(r) - would slow down the particle within a 
characteristic time r. 

Since the physical vacuum does not represent an embed- 
ding medium whose stochastic forces can cause a par- 
ticle to slow down completely, we modify the character 
of the stochastic forces by assuming that they change 
periodically after a laps of « r sec from down-slowing 
"Brownian" to motion enhancing "anti-Brownian" and 
vice versa. The "anti-Brownian" forces act as if the kine- 
matic viscosity would have a negative sign. Hence, the 
corresponding equation of motion has the form 



However, from Einstein's theory of Brownian motion the 
kinematic viscosity (or "diffusion constant" ) emerges as 

k T 

v = — — - (Einstein: ~Ax~Ax] = 2 v At ; (22) 



m 



i,3 = 1,2,3) 



where mo is the mass of the particle under study, 
Axi , Axj are displacements of its position and T is the 
effective temperature of the embedding medium. This 
temperature enters into the derivation as the width of 
the distribution of the random (Gaussian) forces that act 
on the particle apart from the directional down-slowing 
force. Because of the latter there is a down-slowing mo- 
tion that we have already alluded to. The associated time 
constant is denoted by r. Equating the down-slowing 
forces to zero amounts to r — > oo which would yield in- 
finite kinematic viscosity. Hence, there is no alternative 
to our approach. 

Obviously, the physical dimension of the numerator of 
the above fraction in Eq. (|2"2"|) is that of an action, i.e. 
energyxtime. As v appears via V s toch{ r -,t) m Eq. (|2"Tj) 
which is constructed to describe dissipationless motion 
in a "stochastic vacuum" whose effect on a particle can 
only be associated with a new constant of nature, one is 
justified in equating ksTr with A Ti where h = 2irh is 
Planck's constant. Of course, instead of 1/2 there could 
be any other dimensionless prefactor in front of h, but it 
turns out that the numerical results of all quantum me- 
chanical calculations that follow from Eq. (|2"Tj) are only 
consistent with the above choice. Clearly, that choice 
has to be made only once and for all. 
Having thus calibrated the "vacuum-;/' we rewrite 
Eq.dHJ) in the form 



m Q —v(r, t) = F(r) - VV stoch (r, t) - fi(r, t) . (20) 



In Section IXVI we give an example of an embedding 
medium that acts on a test particle by alternating 
Brownian/anti-Brownian forces. 

If we now additionally assume that the temporal changes 
that occur with all quantities in Eqs. (fT7|) and ((20)) are 
slow on a scale of unit length r - which is the stan- 
dard requirement also in diffusion theory - the motion 
of the ensemble will be governed by the arithmetic mean 
of these equations , that is by 



m —v(r,t) = F(r)-VV st och{r,t). 



(21) 



A more detailed definition of the stochastic forces that 
ensure "conservative diffusion" will be given in Section 
IXIII One might suspect that our subdivision into a Brow- 
nian "B "-ensemble and an anti-Brownian "A" -ensemble 
is unnecessarily clumsy and could be avoided at the out- 
set by assuming vacuum forces that neither possess down- 
slowing components nor counterparts that effect motion 
enhancement, but rather consist of random (Gaussian) 
forces whose components form a normal distribution. 



mo—v{r, t) = F(r) - W QP (r, t) 



(23) 



where we have substituted the subscript of V s toch(f, t) by 
"QP" 



V QP = 47^ 



1 (YE 

2 ^ p 



V 2 p 



"quantum potential"(24) 



and wc have set 



2mo 



k B Tr 
m 



(25) 



The "quantum potential" has first been introduced by 
de Broglie^S and later been taken up again by David 
Bohm 28 . Obviously Eq. (f2"3"|) may be viewed as a mod- 
ification of Newton's second law. 

The assumption made above, viz. that all changes of the 
ensemble properties have to be sufficiently slow on a time 
scale of unit length r may raise questions about the va- 
lidity of such a constraint. Eqs. (|23|) and ([24]) will prove 
equivalent to the time-dependent Schrodinger equation 
whose validity is unquestioned at the non-relativistic 
level. Hence, r is obviously sufficiently small within the 
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experimentally tested range of the Schrodinger equation. 
Conversely, as one may conclude then from Eq. (j2"5)l the 
"effective temperature" of the vacuum must be very high 
compared to those temperatures commonly considered in 
applied thermodynamics and astrophysics. 

Fundamentally different from our approach Bohm2£ de- 
rives Eqs. (j2"5)) and (EM)) by choosing the opposite direction 
starting from the time-dependent Schrodinger equation 
which he just considers given. Hence, he does not offer 
any new insight into what makes the motion of a micro- 
scopic particle different from what classical mechanics 
predicts. In the context of Bohm's mechanics Eq. (j23|) 
is frequently cast such that it resembles the Hamilton- 
Jacobi equation. To this end one sets 

v(r,t) = — VS(r,t) and p(r,t) = R 2 {r,t) 

which implies, again for no obvious reason, that v(r,t) 
is curl-free. 

Eq. ([23l) in conjunction with (f24|) then attains the form 



^V 

to 



dS (VS) 2 . , h AR 

1- — + V(r) 

dt 2m y ' 2 m R 



= 0. 



This is equivalent to 



dt 2m [ ' 2to R ' 

and becomes identical with the Hamilton- Jacobi equation 
in the limit Ti — ► 0. However, the connection to classi- 
cal mechanics is far more evident from Eq. ([23|) . which 
reduces to Newton's second law 

d 
at 

as Ti tends to zero. In addition, Eq. l[23|) lends itself to 
a thought-experiment that is particularly illustrative of 
the quantum character of particle motion. 
One starts with setting T% = and assumes that all parti- 
cles of the ensemble commence their motion under iden- 
tical initial conditions. Their positions and trajectories 
will coincide then at any later time. One now lets fi take 
on a finite value. As a consequence of the now occurring 
stochastic forces whose action on some particle is statis- 
tically independent from that on any other particle, the 
particle positions start diverging and form a cloud around 
the formerly common position along the trajectory. The 
particles of the ensemble now reach positions that are 
not accessible under energy conservation. It is hence ob- 
vious that the vacuum provides an embedding medium 
of a "universal noise" consisting of energy fluctuations 
which cause shifts of the individual particle trajectories 
such that the classical momentum and the energy are 
conserved on the average. This is reflected in the expec- 
tation value of the "vacuum force" Fqp = — VVQp(r,i) 
which equals zero: 



J p(r,t)F QP (r,t) d : 



r = 0. 



(26) 



We shift the proof of this equation to Section IXIIII 
Eq. ([26"|) may be interpreted in the sense that the par- 
ticles undergo only reversible scatterings. Figuratively 
speaking, the vacuum keeps track of the energy balance 
and remembers at later positions of a particle departures 
from its classical momentum and energy that occurred 
at previous positions. The undulatory properties of the 
probability density reside in this memory effect which 
gives rise to an unfamiliar non-locality. Hence, from our 
point of view it is illegitimate to correlate these proper- 
ties with a wave-like character of the particle. We 
definitely side with Nevill Mott (1964) who argues: 

"Students should not be taught to doubt that electrons, 
protons and the like are particles. . . . The waves cannot be 
observed in any way than by observing particles. " 



V. THE TIME-INDEPENDENT SCHRODINGER 
EQUATION 

As a first application we discuss the stationary state 
of a particle that is bound to a potential without sym- 
metry elements. Hence, the real-space dependence of the 
potential does not display any distinct direction. That 
means, when a particle of the ensemble appears with a 
velocity Vi(t) in the elementary volume A 3 r around r 
there will always be another particle in that volume with 
approximately the opposite velocity, so that 



v(r,t) 



1 



n(r, t) 



n(r, t) 

Hence, Eq. (f2"5|) reduces to 



(27) 



A mo 



1_ 2 1 /Vp 

P 2 V P 



This is equivalent to: 



4 to 



Vp 
P 



V(r) =0. 



+ V(r)=E, (28) 



where E denotes a constant. Eq. (|28|) represents a non- 
linear partial differential equation in p(r). 
On substituting p(r) by a function ip( r ) defined through 



p(r) = V> 2 (r) 
one obtains because of 



(29) 



Vp = 2ipVip 



1 /Vp 



2 V P 



and 



V 2 p = 2^V 2 V + 2(VV0 5 



--V 2 p = -2 



vV /v 2 f 2 
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a linear differential equation 



If one now makes use of Eq. (I33|) instead of Eq. J29] 



2 mo 



V 2 i) 



V(r) = E that is 



2 m 



X7 2 ip + V(r) V> = E i}) 



(30) 



which constitutes the time-independent Schrodinger 
equation. 



VI. INCLUDING CURRENTS 

For the familiar problem of a particle in a box Eq. ((30|) 
reduces in the one-dimensional case to 



— + k 2 

dx 2 



tp{x) = 



(31) 



where we have set 
A: 

and 



2 2m jp . 



E 



Too 2 
2 17 



F(x) = 





DC 



for < a; < a 
else 



The solutions 



ip(x) 



sin k n x where k n 



/a/2 

may be recast as 

iK*) = 75 hMz)+V>-(z)] 

where 



n; n = 1,2,3..(32) 



■0± {x) 



<p(x) 



k n x 2 • 



In the spirit of our approach the two independent solu- 
tions to the differential equation (|3"Tj) , i/j± (x) , refer to the 
particle moving at a velocity u„ = either to the right 
or (after reflection at x = a) to the left where it is re- 
flected again at x = 0. 

We are thus led to surmise that we have in the general 
case of a freely moving particle 

Mr) = \tp(r)\ e^M and v(r) = — Vip(r) . (33) 

m 

The validity of this conjecture will be shown in Section 

Em 

In a stationary state of the one-particle system in which 
■gjV = but v(r) / we have according to Eq. (fT9|) 
4?v = v-X7v = ^Vi; 2 where we have exploited in advance 
that, according to Eq. ([33"]) . v(r) is curl-free. Hence, in 
the presence of a stationary current Eq. (|28|) contains the 
kinetic energy ^ v 2 as an additional term, that is 



4m 



P 2 



1 fVp 



P 



V(r) 



^V=£.(34) 



P {r) = |^(r)| 2 = U(r)e 



-vp(r, 



(35) 



and substitutes v(r) by V(f(r) the bracketed term in 
Eq. ([34|) becomes 



1 fVp 



4m 



P 



+ i 



h 2 l 9 ft 
5 — tV 2 ^ + 

2 m ip 

■ h 2 




2 771 Q 



vV 



2 rn V V 



Invoking the equation of continuity in the form 

V ■ j = V ■ p v = pV ■ v + v ■ Vp = 

it can readily be shown that the term £[...] on the right- 
hand side equals -mot) 2 . Hence we have from Eq. (f3~4")) 



h 2 l 



2 mo V' 



VV + V"(r) = £, 



that is 



2 mo 



(36) 



as before without a current. 

It should be noticed that ip may well be time-dependent 
even when Vy? is not, that is, we have in general 

¥>(r,t)=¥>o(r) + /(t) 

where f(t) is a real- valued function. In this case the wave 
function ip(r,t) attains the form 

V>(r,i) = V-(r)e l/(t) where $(r) = |^(r)| e 4Vo(r) (37) 
and hence, its time-derivative may be cast as 
d 



ih^(r,t) = -hfiP(r,t). 



(38) 



Since f(t) is primarily unspecified and —hf possesses the 
dimension of an energy the latter may justifiably be iden- 
tified with the energy E which is the only energy-related 
constant characterizing the wave function of the system: 



hf = E; that is if(t)=--Et. (39) 

h 



As a result, we have from Ea. ([3T|) 

ip(r,t) = ^(r)e~^ Et 



(40) 
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for a wave function in a stationary state. Furthermore, 
we have from Eqs.([56]). ([55]) and (03) 

- ^— V 2 ^(r, t) + V(r) 1>(r, t) = %h ^(r, t) (41) 

which constitutes the time-dependent Schrodinger equa- 
tion. Its validity is here still restricted to stationary sys- 
tems, but it will be shown in Section|X]that it retains this 
form also for non-stationary systems. However, in order 
to achieve this consistency, one has to introduce the neg- 
ative sign in Eq. (|39[) which seems to lack reason and can 
actually not be justified without reference to Section IXl 

VII. THE VELOCITY POTENTIAL AND 
PHASE UNIQUENESS 

We rewrite Eq. (|2"3")) in the form 

m -v(r,t) = -VP{r,t) (42) 
at 

where 

P(r,t) = ^-[V(r) + V PQ (r,t)], 

and we have made use of Eq. (Til)|) defining the "hydrody- 
namic" or convective acceleration 

± v (r,t) = -^v+(v.V)v. 

In hydrodynamics Eq. (j4"2")) corresponds to the Euler equa- 
tion of perfect (frictionless) fluids and constitutes the 
starting point of Helmholtz's theory of vortices. Thom- 
son's more elaborate analysis on vortices?- builds on 
Helmholtz's considerations. We confine ourselves here 
to reporting only the general ideas as far as they directly 
concern the present theory. 

If we set uj — V x v for the curl of the ensemble average 
of the particle velocity, we have from Ea. (TTB| 

v 2 

(v ■ V) v = V v x uj . 

v ; 2 

We now form the curl of Eq. (j42|) and use this expression 
together with Eq. (fT9|) . The result may be cast as 

-| uj(r, t)-Vx [v(r, t) x tD(r, t)] = , (43) 

where we have used V x VP = and V x Vi> 2 = 0. 
One recognizes from Eq. (|4"3")) that Jj cJ(r, t)\ t =o becomes 
zero for some chosen time, which we here equate to zero 
for convenience, if uj(r,t)\t=o — at that time. Forming 
the time derivative of Eq. (f43|) and setting again t = we 
see that the second time derivative of uj(r,t) vanishes as 
well. This can be carried further to any higher order of 
the time derivative. Hence, the system stays curl-free if 



it is curl-free at t = 0. We now consider an ensemble of 
free particles (F(r) = 0) when U — 0. They may start 
their motion at t = at the same point in real-space and 
with the same momentum po = rriQVQ. If one allows Ti 
to attain its natural value, the particle positions diverge 
and form a point cloud. Outside this cloud there are 
no particles and therefore v(r,t) = 0. Since the ensem- 
ble does not exchange momentum with the vacuum on 
the average and consequently no angular momentum, we 
have everywhere within the space of normalization 

V x v(r,t) = uj(r,t) = Vr,i. (44) 

If one now turns on some (physically realistic) potential 
V(r), weighting it with a smooth switch function from 
zero to one, starting at t = to, the velocity distribution 
v(r,t) for t > to will now change differently, of course, 
but because of Eqs. (f4"3"]) and (T4"4"]) for t = to, we have as 
before uj(r,t ) = and ^ w(r, t)| t=to = which again 
applies to any higher order time-derivative at t = to. 
We thus arrive at the conclusion that an ensemble whose 
equation of motion is given by Eq. (|42p is curl- free. In 
other words, v(r,t) possesses a potential ip(r,t) which 
we express in the form 

v(r,i) = — V(/3(r,i). (45) 
too 

Because of the prefactor Ti/mo the function <p(r,t) be- 
comes dimensionless. Eq. (l45l) may equivalently be cast 
as 

uj{r) = / v(r') ■ dr' (46) 

" Jr 

where we have omitted the time-dependence in confining 
ourselves to a stationary state situation. As in the theory 
of perfect fluids there may be singular vortex lines which 
occur if V(r) possesses axial or spherical symmetry. A 
vortex line then defines an axis of quantization. The 
latter may be regarded as the boundary line of a semi- 
plane. Even in the presence of a vortex line, can (p(r) be 
defined such that it remains unique if one only stipulates 
that the starting point of the line integral in Eq. (|46)) . ro, 
lies on one side of this semi-plane and that the path along 
which the integral is performed never crosses that semi- 
plane. The point ro may be chosen at will. In general, 
(p(r) will now be discontinuous at the semi-plane. The 
ensuing section deals with this particular problem. 

VIII. QUANTIZATION OF ANGULAR 
MOMENTUM 

The primary objective of this section is to disprove 
Wallstrom's notable objectionii against Madelung's con- 
viction, also held by other theorists of this school of 
thought, that Newton's modified second law (|2"3")) is equiv- 
alent to the time-dependent Schrodinger equation which 
we shall derive below. In so doing we have to exploit the 
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uniqueness of the velocity potential shown in the preced- 
ing section. By contrast, in standard quantum mechanics 
the time-dependent Schrodinger equation is regarded as 
given. It is customarily converted into the equation of 
continuity 

{f Vf/--^Vf }] = 



L 1% niQ 



to show that the bracketed expression has to be inter- 
preted as the current density j(r,t). This conclusion is 
only legitimate if j(r,t) has been proven to be curl-free 
which, however, is only tacitly presupposed. Inserting 



(47) 



^(r,t) = \i;(r,t)\e^^ 
into the bracketed expression yields 



j(r,t) = \^(r,t)\ 2 _Vv(r,t), 

— v{r.t) 

as a consequence of which one obtains Eq. (|45|) . If one is 
dealing with a stationary state whose velocity field con- 
tains a vortex line, e. g. an excited state of a hydrogen 
electron possessing an orbital momentum, we have 



<j> v(r) ■ dr =/= 



(48) 



for any path encircling the vortex line (=quantization 



axis). On inserting here v — — one obtains 

/ ° mo ' 



V<p(r) ■ dr — ip(r) — ip(r ) ^ 



(49) 



where r and r^ are two points facing each other across the 
semi-plane, introduced in Section [VIH at an infinitesimal 
distance. Thus, in general the phase of the wave func- 
tion, and consequently the wave function itself, will be 
discontinuous at the semi-plane as opposed to p(r) and 
j(r) which may be presupposed to be smooth functions 
everywhere. 

Clearly, as follows from Eq. (|4"T|) , ip(r) remains continuous 
at the semi-plane if 



tp(r) — <^(j"o) — 2to7t where 



integer . (50) 



But there is no immediately obvious reason why one 
should require tp(r) to be continuous because only p(r) 
and j(r) can be regarded as reflecting physical properties 
of the system. We are hence led to conclude that with- 
out an additional argument neither our derivation nor 
standard quantum mechanics yields a justification of the 
proven relation 

itlq <j> v(r) ■ dr = 2mn h — mh where m = integer (51) 

which comprises Eqs. (|4"5"|) , (|4"5|) to (|5D|) . This has already 
been pointed out more than 75 years ago by Paul: 26 



and Born and Jordan^!. As opposed to these consid- 
erations Wallstrom states in his paper-ii: "To the best of 
my knowledge, this condition fEq. tfSTj) ) has not yet found 
any convincing explanation outside the context of the 
Schrodinger equation". 

This is definitely incorrect: within that context the as- 
sumption of continuity fEq. ([5(1)0 has to be justified by 
an additional argument as well. 

What else necessitates then the continuity of ip(r) every- 
where? 

We consider two states, V> mi (r) and -0 m2 (r), which are 
solutions to the time-independent Schrodinger equation 
for a spherically symmetric potential. The associated en- 
ergies may be denoted by E mi and E m2 , and the spherical 
coordinates by r, 9, a with 9 = 0, it defining the quantiza- 
tion axis (=vortex line). We assume that there is a weak 
magnetic field along this axis to lift a possible degeneracy. 
The solutions have the general form 



1 



2tt 



Fi 



ki k 2 k 3 



(r,6)e 



i k$ a 



where k\,k,2, ks are real numbers and 

tp(r) = ip(a) = fc 3 a = m 1/2 a . 



(52) 



Moreover, F kl k2 m 1/2 (t, 9) denotes real- valued functions 
whose square is normalized to unity. Specifically we dis- 
tinguish two solutions which differ in their /^-values and 
solve the time- dependent Schrodinger equation 

i, mi/a (r,t) = -}= F kl k2 rni/2 (r, 9) e i( m ^ 4) 
v ^ 

where 

^ m l/2 — Fra-^j 2 lh . 

Since the time-dependent Schrodinger equation is linear, 
the function 

4>[r, t) = ~h [crm F kl k 2 mi (r, 0) ^ ™ «> 



~\~Cm 2 F kl k2 rrL2 (r, 9*) 6 



,( m 2 a—CjJmn i)l 



constitutes a solution as well. The constants c mi/2 may 
be assumed real-valued. We now form 



^)l — c m 2 Cm 2 I-rni '< 

=p(r,t) 



1 

2^ 



2tt 



; i(m 2 -mi)a ^ gi(w mi -w m2 ) t _|_ 



where 



I-mi m 2 = / F kl k2 mi (r, 9) F kl k2 m2 (r, 9) r dr sin 9 d9 . 



According to Eq.Q the real-space integral of p(r, t) must 
be unity at any time. This is obviously only possible if 
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TO2 — Tii is an integer. Since mi may attain the particu- 
lar value zero, corresponding to a state without circular 
current, it follows then that to 2 must be an integer itself 
and hence, if one invokes Eq. ([5^ 



too j> v(r) ■ dr = h <j> V(p{r) ■ dr 

/•27T Q 

= h — — (to2 a) da = TO2 h 
Jo 9a 

where to 2 is now proven to be an integer number, in 
agreement with the conjecture (fSTj) . 
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FIG. 1: Interference of trajectories 



IX. AN INSTRUCTIVE OBJECTION AND 
QUANTUM BEATS 

An apparently serious objection against a stochas- 
tic foundation of quantum mechanics along the lines 
of the preceding sections goes back to Mielnik and 
Tengstrand 31 . The authors refer to an experimental 
setup as sketched in Figure [1] where the test particle en- 
ters from a distant source on the left-hand side and is 
kept within a tube that extends up to a screen on the 
right. The tube contains an impermeable partition that 
completely seals off the upper part (A) from the lower 
part (B). It possesses a limited, but macroscopic length 
of, say, 10 cm. The authors argue that according to con- 
ventional quantum mechanics the incoming wave would 
split up into an upper and totally independent lower por- 
tion. Yet both portions retain their capability of inter- 
fering with each other when they merge again within the 
area C and beyond. However, if the wave portions are re- 
placed by the set of irregular trajectories which stochas- 
tic quantum mechanics claims to be an equivalent of, 
it seems to be very unlikely that stochastic-force con- 
trolled trajectories can preserve information over so long 
a distance as well as waves. This criticism amounts to 
perceiving the preceding derivation of the Schrodinger 
equation from Eq. (|23p as ill-founded or even erroneous. 
It is just the solution to the Schrodinger equation for 
the particular setup around which the authors' consid- 
eration revolve. On the other hand, it is easy to ver- 
ify the validity of the derivation. There is simply no 
step where one may be in doubt. But one has to keep 
in mind that the solutions i/j(r) = \ip{r)\e lv ^ to the 
Schrodinger equation provide only information on ensem- 
ble properties and not on a particular trajectory that is 
a member of the ensemble under study. For example, 
the velocity v(r) = ^ V<p(r) at some point in the area 
marked C represents such an average over all trajectories 
of the ensemble running through that point. This en- 
semble defines the probability in which direction a par- 
ticular particle that has arrived at C, e. g. along the 
"A" -trajectory, will move further. (S. Figure [IJ lower 
panel.) This is analogous to considerations we shall dis- 
cuss in the context of the Smoluchowski equation (Sec- 
tion [XXXlT]). The properties of the ensemble are just an 



image of the property of the vacuum fluctuations to en- 
sure the absence of dissipation. This manifests itself in 
the fact that v(r) is curl-free as in ideal fluids. An indi- 
vidual particle that has moved along the "A" -trajectory 
and arrives at C "feels", so to speak, the possibility of 
a "B" -trajectory. As stated above, it continues its tra- 
jectory depending also on the family of "B" -trajectories 
running through C. If the partition in the tube would be 
elongated and the point C correspondingly shifted to the 
right, irrespective of how much, the "A"- and "B" subset 
of trajectories would now be different, but the scattering 
probability at C of a particle that has moved along an 
"A" (or "B ")-trajectory would still be influenced by possi- 
ble "B" (or "A")-trajectories. Furthermore, if one would 
place some electrostatic array into the upper part of the 
setup which would cause a spatially varying electrostatic 
potential, the "A" -trajectories would change accordingly 
and give rise to a different interference pattern within the 
"C" -range. To make the surprising content of this obser- 
vation even more striking we consider a situation where 
one has particles enter the setup one by one from the left 
so that only one particle traverses the setup at a time. 
First, we switch the electrostatic array off so that there 
is no extra potential along the "A"-trajectory. If one has 
placed a detector, an electron multiplier, for example, 
at some position r screen on the screen, it would monitor 
the incoming electrons at a certain rate. These electrons 
come either along an A"- or a "B "-trajectory. Once the 
extra potential has been turned on, the count rate at 
f screen changes. Although an electron may have moved 
along the unmodified "B "-portion of the setup, it feels 
the modification of the "A" -portion when it arrives in 
"C" . As explained above, this is due to the change of the 
vacuum scattering probability at C. Electrons that have 
arrived at some elementary volume within C and have 
so far preferentially been scattered into r screen are now 
also scattered to other positions on the screen, thereby 
changing the count rate at r screen . 

If the electrostatic array in the "A" -portion would sim- 
ply consist of two planar parallel grids perpendicular to 
the average particle motion, and if one applies an accel- 
erating voltage V between the grids, the particles' kinetic 
energy eo increases by an amount Ae — eV where e de- 
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notes the particle charge. The wave function tp(r screen , t) 
at the screen is the sum of the "A"- and "B" -related con- 
tributions: 



1 

V2 



(Vscreen) e~™ A 1 + Vs^creen) e~^ B 1 (53) 



^ A, B {r screen) = ^e' kA ^ r - 



% u A = e + Ae ; Ti ui B = e , 



where 



and 



and with 1/y/N denoting an appropriate normalization 
factor. The function ip(r,t) solves the time-dependent 
Schrodinger equation (|4ip for the particular array under 
study. If we introduce e = eo + ^ Ae we may cast h kj^iB 
as hk A / B w V 2 TOoe (l ± \ ^r) if 4^ << 1 where to 
denotes the rest mass of the particle. Eq.([53| can then 
be rewritten 



V2 



ipA/B(r scr 



i\ — _L_ p'C" 



-i(Ak r sc 



It) 



-Auj t) 



where % k = y/2mo e , Afc = k A — ks and ?i Aw = i Ae. 
Hence we have for the current density j(r,t) oc count 

rate at fscreen 



3{r s 



>*) = W wIV'A/Bt'Wen,*)! 2 



N m 



mo 

[1 + cos(2Afcr s 



-Ae*)] 



That means: the count rate oscillates at a period of 
T = . This most surprising effect of "quantum beats" 
has clearly been observed by Rauch and collaborators (s. 
Badurek et alt^S) who used spin polarized neutrons in- 
stead of electrons. The energy change Ae in the "A"- 
section of the setup was in that case imparted to the re- 
spective neutron by flipping its spin within a spatially 
confined magnetic field along the "A" -trajectory. (In 
practice one used a spin flipper also in the "B" -portion 
of the setup where the corresponding magnetic field was 
slightly lower than in the "A" -portion so that Ae referred 
to the difference of two spin flip energies in this case.) 



X. THE TIME-DEPENDENT SCHRODINGER 
EQUATION 

In the most general case v and hence tp are time- 
dependent. As already pointed out in Section I VII the 
substitution of p(r) has to be modified then in the form 



(54) 



which was introduced by Madelung in 1926^i The ±- 
sign requires a comment. As discussed in Section 2, 



p(r,t) will generally be presupposed as a smooth func- 
tion. The zeros of p(r,t) pose a particular problem that 
occurred already in Section 3, but was not explicitly men- 
tioned. The admissible type of zeros limits the set of 
functions p(r) that can be mapped onto ip{ r ) according 
to Eq. (f35|) . For simplicity we confine ourselves to the 
time-independent case and assume that the zeros of p(r) 
lie on the faces of a rectangular parallelepiped defined by 
the equations x v = x^q with x v and v — 1,2,3 denot- 
ing Cartesian coordinates. Hence close to x v = x v q and 
perpendicular to the respective face the density varies as 
(x„ — Xvo) 2 . Since we have everywhere p(r) > its square 
root varies as \x v — x v $\ and thus would not be differen- 
tiable at x v — x v q. In defining the map p{r) — > ip(r) 
one is forced hence to choose the positive sign in front of 
-J p(r) outside the rectangular parallelepiped if one has 
chosen the minus sign inside (or vice versa) to ensure that 
ip(r) stays differentiable across the face of the rectangu- 
lar parallelepiped. Hence, mapping functions p(r) onto 
differentiable functions ip(r) is only possible if the zeros 
of p(r) subdivide the space of volume V into cells without 
leaving empty space. At first sight it appears that this 
limitation in the set of admissible functions p(r) consti- 
tutes a serious drawback of the entire concept. One has 
to bear in mind, however, that the functions ip(r) are 
not determined as a map of p(r) but rather by solving 
the Schrodinger equation ([3H)l which has been the objec- 
tive of the derivation. Physical meaningful solutions to 
Eq. (|30|) have automatically the required spatial structure 
of their zeros. 

We now move on to derive the time-dependent 
Schrodinger equation under the supposition that the 
above considerations apply to the time-dependent case 
as well. 

If one uses instead of Eq. (|2"3")) the arithmetic mean of the 
original Eq.([7]) and its "anti-Brownian" analogue where 
the sign of v and u{r,t) is reversed, one obtains 

d ft 1 

— v + (v ■ V) v - (u ■ V) u H Au= — F(r) (55) 

at 2 rrio m o 

This can be simplified in the form: 







1 



1 



1 



-v = W - - Vtr + - Vu A - 



dt too ' 2 2 2 Too 

where we have made use of the relations 



Am, (56) 



1 2 1 2 T 1 

v ■ Vv — -Vt> : u ■ Vm — -Vm and v — . 

2 2 2to 

With the first two equations it has been observed that v 
and u are curl-free. On differentiating u with respect to 
time and using Eq.® one obtains 



d 



Ti „ ,d p , . 
V(a7//») 



dt 2 too d t 

Invoking the equation of continuity 

dp/dt + V ■ (pv) = 



(57) 



(58) 
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that is 



dp/dt + pW -v + v Vp= 



4t//9 can be replaced by— V ■ v — v ■ - Vp which yields 



U _ ,dp , . h _ ,_ . 
7; V (-^- /p = V (V • v) 

2m y dt /H ' 2m v ; 



h 1 



2 mo p 



(59) 



Using Eq.© we may substitute the expression in the [...] 
brackets on the right-hand side by v ■ u. Hence Eq. ([59 
takes the form 



d_ 

dt 



2 mo 



V(V • v) - V(w ■ v) . 



(60) 



On multiplying the equation of motion (|56p by the imag- 
inary unit i and subtracting Eq. (|60p we obtain 



d_ 

di 



( — u + iv) 



Vu 2 



l lot 

WV--Wv 2 

m 2 2 

i Am - V(V ■ v) + V(tt • v) 



2 mo 



2 mo 



which after reordering the terms on the right-hand side 
becomes 



■Vf-«- 



d r 

+7^-V[V- (-u + iv)] - — VV. (61) 
2 mo mo 

Here we insert Eqs. fTB")) . (l4"5")) and in the form 

+ = — Vln^/Vpo] ■ (62) 
m 

After interchanging the operators 9/9 t and V one ob- 
tains 



V 



V 



n i dip 

mo ip dt 



V [jVip)\-—V 

,W / too 



If the gradient of some function equals that of another 
function the two functions can only differ by a real-space 
independent function of time which we denote by f3(t). 
Hence, if one divides the above equation by the imaginary 
unit the result may be cast as 



h l dip i h 2 



1 mo ip dt 



2 ml 



V-ip(t) (63) 

m 



One can now make use of the identity 



and multiply Eq. (|63p by —rriQip. This yields 
ft 2 V 2 



9^ 



2 mo 



V> + V ip + j(t) ip 



(64) 



where 

7(f) = im Q (3(t) . 
If ip(r,t) is replaced by ip(r,t) defined through 



ip(r, t) — ip(r, t) exp 



~ / iV) at 1 



to 



Eq. ([M|) becomes an equation for ip(r,t): 



diHryt) 
dt 



2 mo 



V(r) 



iP(r,t). 



where 



p = —ifiV . 



(65) 



(66) 



The two functions ip(r,t) and ip(r,t) differ only by a 
time-dependent phase factor without physical relevance. 
Only the functions 



p{r,t)=iP*(r,t)iP(r,t) 
and the current density: 

j(r,t) = p(r,t) — V<p(r,t) 
m 



(density) (67) 



(68) 



refer to relevant quantities of the system which obviously 
do not depend on this phase factor. For this reason we 
may set ^(t) = 0, that is replace ip(r,t) in Eq. (f6"5)) by 
ip(r,t) without loss of generality. To simplify the nota- 
tion we introduce the so-called Hamiltonian defined by 



H 



P 

2 mo 



V(r). 



(69) 



Eq. ([65|) then takes the familiar form of the Schrodinger 
equation 



in^m=H(r)iP(r,t). 



(70) 



The first order time derivative on the left-hand side can 
be tracked down to the acceleration (d/dt) v in Newton's 
modified second law (|55|) . 
Using 

iP(r,t) = \iP(r,t)\e lvir ' t} 

and inserting this into Eqs. ([6T|) and (f6"8"|) one obtains the 
familiar expression 

j(r,t) = p{r,t)v(r,t) = 

[iP*(r,t)ViP{r,t)-iP(r,t)ViP*{r,t)] (71) 



2i mo 
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which on real-space integration and multiplication by too 
yields 



(«(<)) = / r (r, t) p ^(r, i) d J r = (p> (72) 



to • 



where t/)(r,t) has been required to satisfy the usual 
boundary conditions at the surface of the normalization 
volume. Because of Eq. (|72"|) one is justified in terming p 
"momentum operator" . 

In Bohm's version of quantum mechanics 28 Eq. (|71|) is re- 
cast to define the velocity field 



v(r, t) = — 



/ V^(r,t) 
m ~ V 



The streamlines of this field are interpreted as true par- 
ticle trajectories. From our point of view this appears to 
be rather absurd because the explicit r-dependence of v 
comes about by forming the ensemble average over the 
(in principle infinite) family of true trajectories as defined 
in Eq.Q. Bohm's definition of v as describing the true 
velocity of the particle leads inescapably to strange re- 
sults, notably with stationary real- valued wave functions 
ip(r) for which v(r) = 0. Hence, the particle appears to 
be at rest although the kinetic energy of the particle 



(T) 



P 



tp*(r) V , ( r ) d r 

2 toq 



{Pi 
2 toq 



(73) 



is definitely different from zero. 

The time-dependent Schrodinger equation represents the 
center of non-relativistic quantum mechanics. Funda- 
mentally different from the present approach where it 
is derived from a new vacuum concept, in conventional 
quantum mechanics it falls out of the blue, and this ap- 
plies to Bohm's theory as well. As the latter associates 
the pattern of smooth stream lines with the set of true 
particle trajectories, it is forced to explain the probabilis- 
tic character of the information contained in ip(r, t) by 
an additional "quantum equilibrium"- hypothesis. It is 
therefore hard to see that anything can be gained by "go- 
ing Bohmian" . The "process of measurement" in which 
a particle moves from a source to the detector where it 
fires a counter, is in our view described by one of the ir- 
regular trajectories which is terminated at the detector. 
Due to the stochastic forces that cause this irregularity, 
the information on the ensemble properties is naturally 
probabilistic. 

A frequently raised objection against Bohm's theory con- 
cerns the asymmetric way in which it deals with the par- 
ticle's real-space position and its momentum. In fact, 
the real-space position r plays a pivotal role in Bohm's 
theory compared to the other observables which are "con- 
textualized" by resorting to the wave function tp(r,t) 
that solves the Schrodinger equation for the system un- 
der study. By contrast, in our approach the ensemble's 
i-th particle position Vi and its velocity V{(t) enter into 
the theory as autonomous quantities. This is reflected 
in the occurrence of two independent functions p(r,t) 



and v(r,t) 



■V<p(rt). It is this pair of informa- 



tion p(r, t) ; ip(r, t) that necessitates the description of 
the one-particle system by a complex-valued function 



XI. THE UNCERTAINTY RELATION AND THE 
ISSUE OF "MEASUREMENT" 

By performing a Fourier transform on ip(r,t) 



C{k,t)e lkr d 3 k (74) 



(2tt) 3 / 2 

E cis. ([72]) and (|73|) may alternatively be written 



where 



(P) = I ^*(r,t)pt/j(r,t)d 3 r = 
J C*(k,t) (hk)C(k,t) d 3 k 
(P 2 ) = I r(r,t)p 2 ^(r,t)d 3 r = 

C*(k,t)(hk) 2 C{k 7 t)d 3 k (75) 



(76) 



C*(k, t) C(k, t) A 3 k = P(k, t) A 3 k 



describes the probability of the particle possessing a mo- 
mentum that lies within A 3 fc around k in the fc-space. 
We temporarily label the coordinate-components of the 
particle in the two spaces by an index v\ v — 1,2,3. 
The mean square departures of the position coordinates 
x v and k Vl respectively, from their arithmetic means x v 
and k v are given by 



{{x v -x u f) t = I i}*(r,t)(x v -x l ,f^{r,t)d 3 r 



and 



((** -h?)t = / C*(k,t)(k„-k„) 2 C(k,t)d 3 k. 



Since C(k,t) is the Fourier transform of ip(r,t) we have 
as a fundamental mathematical theorem 



{{x v -x u ) 2 ) t {{k v -k v f) t > 



1 



that is 



Following the standard notation by setting Ax,, 



(77) 



a/ {{x v - x v ) 2 )t and Ap v = \/ {(h k v -% k u ) 2 ) t = 
Ax u A Pv > - (78) 



y/ {(p— (p)) 2 )t the latter relation may be cast as 

ft 
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which constitutes the celebrated uncertainty relation. It 
is commonplace to interpret this relation, loosely speak- 
ing, by saying: "momentum and position of a particle 
cannot be measured simultaneously with any desirable 
precision" . 

From our point of view it does in no ways refer to any 
measurement on the position or momentum of the par- 
ticle in question. It is nothing more than the theorem 
Eq. (|77|) on the product of two quantities that are inter- 
connected by a Fourier transform. Furthermore, since 
this relation is - besides the Schrodinger equation - just 
another consequence of our concept, it cannot possibly 
conflict with the existence of trajectories which consti- 
tute a fundamental element of that concept. 
Eq. ([78|) is considered ground-laying for the Copenhagen 
interpretation of quantum mechanics. The latter is based 
on the conviction that it is the measurement that causes 
the indeterminacy in quantum mechanics and necessi- 
tates a probabilistic description of microscopic mechan- 
ical systems. In a highly respected article^ Heisenberg 
gives a revealing example of such a measurement. To 
pinpoint an electron moving along the x-axis within an 
experimental setup he considers a 7-ray source, that il- 
luminates the electron beam, and a hypothetical 7-ray 
microscope that possesses a sufficiently high resolution 
in detecting the position of that electron up to an error 
of Ax. He demonstrates that the 7-ray photon that "hits 
the electron" and is subsequently scattered into the mi- 
croscope, transfers a momentum Ap x to the electron so 
that 

Ax Ap x w h . (79) 

The above result reflects only a property of the micro- 
scope 

Ax Ak x sa 2 tt 

which interrelates the resolved linear dimensions 
Air = A/ sin a of an object and the admissible maximum 
angle a required to ensure that the scattered wave (of 
wavelength A) is still captured by the front lens of the 
microscope, and Ak x = k sin a which describes the 
A^-change of the wave vector of the scattered wave. 
But this interrelation expresses only the content of 
Eq. (|77[) in a different form. The measurement, how- 
ever, is completely fictional for two reasons. Firstly, 
imaging systems within that regime of wavelength are 
for fundamental reasons unfeasible. Secondly, different 
from the picture insinuated by Heisenberg's phrasing, 
the interaction does not take place as an instantaneous 
collision process where a point-like particle (the photon) 
hits another point-like particle, the electron. Instead 
the transition probability of the electron for attaining 
a different momentum is given by the mod squared of 
the transition matrix element M opt , a real-space integral 
that extends over a range of many light wave lengths in 
diameter. Moreover, the transition is not instantaneous 
but rather takes some time of the order fi/\M opt \. Within 



this transition time the electron covers a distance Ax' 
which has nothing to do with Ax in Eq. (|T9")) . Other 
examples of "measurement" , e. g. diffraction at slits 
of a certain width Ax show even more directly that 
the probabilistic information on the (non-relativistic) 
motion of a particle is exhaustively described by the 
Schrodinger equation and boundary conditions for tp(r), 
and hence this information merely reflects our vacuum 
concept, irrespective of whether or not results on the 
diffraction are verified by measurements. 
The host of considerations invoking the uncertainty 
relation (|75|) refer to situations where a particle is 
located within an interval Ax and one interprets this 
confinement of the particle indiscriminately in terms 
of a "measurement" of its coordinate x with limited 
accuracy. One concludes then from the uncertainty 
relation that Ax correlates unavoidably with a variance 

Ap 2 of its momentum such that Ax Ap x « l| where 
Ap x = \J Ap 2 . In reality neither a measurement on Ax 

nor on Ap 2 is truly executable. The uncertainty relation 
merely states that a solution of the one-dimensional 
Schrodinger equation for a particle in a box of length 
Ax yields a ground state energy AE — where 

Ap 2 = (fi-£^) 2 - Hence one obtains simply as a conse- 
quence of solving the Schrodinger equation for that case 
"without observer" (!) Ax Ap — ttTi. One cannot help 
but quote John Bell's question phrased in his stirring 
article "Against Measurement"— 

"What exactly qualifies some physical systems to play 
the role of 'measurer'?" 

The above considerations are in line with a discussion of 
Heisenberg's paper by Wigner— . 



XII. AVERAGING OVER THE TOTAL 
ENSEMBLE 

In forming the arithmetic mean of the two equations 
(fT7| and ([20]) we omitted to mention a problem that we 
wish to discuss here in more detail. 

We temporarily decompose the entire ensemble consid- 
ered so far into a "Brownian" and "anti-Brownian" sub- 
ensemble, each characterized by the associated stochastic 
forces and comprising an equally large number of mem- 
bers. Accordingly we distinguish the velocities v(r,t) 
and the densities p{r,t) in the respective sub-ensembles 
by subscripts B (for "Brownian") and A (for "Anti- 
Brownian"). If the velocities in these two equations agree 
at a time t, they are definitely different at a later time 
t + At. Yet forming the arithmetic mean of the two equa- 
tions can only lead to the same average - which we could 
recast as "Newton's modified second law", Eq. ([2Uf - if 
the two velocities Vb,va and the densities Pb,Pa agree 
also at t + At and any later time. At first sight the lat- 
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ter appears to be irreconcilable with the former. One 
has to recall, however, that our subdivision of the entire 
ensemble into sub-ensembles B and A represents only a 
simplifying model for the actually occurring reversible 
scatterings. In the real system the stochastic forces of 
the J5-type become automatically forces of the ^4-type 
and vice versa within the characteristic time r so that 
the change of the velocity At; in either sub-ensemble is 
[Av A + Av B }/2 within a time span At ^> r which, how- 
ever, must be small compared to time intervals within 
which the quantities of interest change sizeably. The sit- 
uation is similar to that encountered in diffusion theory 
where we have 



dp 
Of 



vAp. 



This equation is obtained from the equation of conti- 
nuity for v B = u and u = —l/Vp/p with the latter 
equation based on similar considerations as the deriva- 
tion of Eq. ([7]) invoking Einstein's law ([22]) which implies 
At » t. The above equation of diffusion hence describes 
changes that are actually defined only on a coarse grain 
time scale and its validity is confined to changes that are 
sufficiently slow on that time scale. As we have already 
discussed in Section HVl this is also the assumption un- 
derlying our derivation of Newton's modified second law 

HMD- 

We temporarily rewrite the two equations (fl7]) and ([20]) 
for an - in that sense - "appropriately long, but suffi- 
ciently short time interval" At in the form 



Av B/A (r,t + At)=R B/A (r,t) At 



where 



R 



and 



B /A(r,t) = -J- (-V[V(r) + Vop(r,t)]±6(r,t)) (.80) 







du 
~dt 



+ (v ■ V) u - (u ■ V) v + v A v . 



Here we have already used Vqp instead of V a toch , but still 
denoted the prefactor of Av by v to demonstrate that f2 
(and consequently u) changes sign when v changes sign. 
It should be noticed that according to Eq. ([24]) Vqp has 
the property 

V QP (p B {r,t)) = V QP (p A (r,t)) = V QP (p{r,t)) since 
at the time t under consideration we have p A (r,t) — 
PB(r,t) = \p{r,t). 

After the elapse of a time At within which each of the 
N/2 particles in the two sub-systems has changed its af- 
filiation (B from A or vice versa) we have 

v B /A(r, t + At) = v B/A (r, t) + Av(r, t + At). 

where 



Av(r,t + At) = AR B/A {r,t + At) = 

— (-W[V(r) + Vp Q (r,t)}) At (81) 
m 



with AR B / A denoting a time average. 

> We consider Eqs.f lHOj) and flHTl) as implicitly 
defining "motion under reversible scattering" <c. 

If we now form the average over the entire ensemble we 
get 

v(r, t + At) = ~v B (r, t + At) + -v A (r, t + At) = 

v(r,t) + — {-V[V{r) + V PQ (r,t)]) At. 

Thus we have 

v B (r,t + At) = v A (r,t + At) = v(r,t + At) . 

We want to demonstrate that the densities behave anal- 
ogously. For this reason we resort to the equation of 
continuity ([58"]) which holds for each sub-ensemble 

dp B / A 



dt 



V • (pb/aVa/b) = 0. 



(82) 



It describes the conservation of the number of particles 
in each of the two subsystems. We conclude from this 
equation that ps(r,t) = p A (r,t), if p B (r,t) = p A (r,t) 
and v B (r, t) = v A (r, t). If one differentiates Eg. ([82]) with 
respect to time and uses v B (r,t) = v A (r,t) = v(r,t) 
as a result of the preceding considerations, we may 
conclude p B (r,t) — p A (r,t). One can carry this con- 
clusion further to any order of time derivative. Thus, 
the Taylor-expansions of p B (r,t + At) and p A (r,t + At) 
agree for any length of the time interval At if p B ,p A 
and v B ,v A agree at time t. 



XIII. CONSERVATIVE DIFFUSION 



We want to prove the validity of Eq. ([26]) which consti- 
tutes a necessary condition for the preservation of classi- 
cal motional behavior on the average. To see this more 
clearly, we first consider one particle (the i-th.) in the cube 
A 3 r around r acted upon by the external force F(r) and 
the stochastic force F si (t). According to Newton's sec- 
ond law we have 



^m ^(t) 



F(n) + F si (t) . 



If we sum this equation over the n(r,t) particles con- 
tained in A 3 r, divide by N and form ensemble averages 
similar to Eqs.|T]) and we obtain 



a i 

di m °N 



E ««(*) 



i=l 



-n(r,t) v(r,t) 



i 

N 



E 

i=l 



F(n) 



=n(r,t) F(r) 



n(rj.) 

- E ^iW ■ ( 83 ) 



=n(r,t) F s (r,t) 
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Here the summation runs over all particles in the cell 
irrespective of whether they belong to the first or second 
sub-ensemble. 

The idea of "conservative diffusion" implies that the N — 
J2r r n ( r >t) particles of the entire ensemble do not feel 
a stochastic force on the average although F s (r,t) does 
locally not vanish in general. Thus, F s (r,t) is required 
to have the property 

^Ht^- F s (r,t) = J p(r,t) F QP (r,t) rf 3 r = 0V^4) 

=F s (r,t) 

yields after summation over 



as a result of which Eq. 
all elementary cells 



d irA n(r,t) , . d 



N 



dt 



p(r, t) mo v(r, t) d 3 r 



=<p(t)) 



E 



n(r, t) 
N 



F(r) 



=j p{r,t) F{r)d 3 r={F) 

We thus obtain as a consequence of the required property 
of F s (r,t) 



We rewrite the integral over the second term on the right- 
hand side using Gauss' theorem 

/ p V ■ (Vu) d 3 r = / V ■ (pVu)d 3 r- Vp-Vud 3 r. 
Jv s « ' Jv Jv 



=Au 



= fj. p Vit-d 2 r 



We assume that p(r,t) differs sizeably from zero only 
within a volume that lies completely within the finite 
space and drops sufficiently fast to zero toward infinity 
so that the surface integral vanishes. Using again Eq.© 
we hence arrive at 

-J Vp-V«d 3 r=^p J p{u-V)ud 3 r 

which shows that, in fact, the right-hand side of Eq. (|86|) 
equals zero. Thus the expectation value of the right- 
hand side of "Newton's modified second law", Eq. ([23|) . 
becomes equal to (F). However, we have on the left- 
hand side (^mo d) instead of 4(mo!)). Nevertheless, 
the two expressions are equal as follows from multiplying 
4? mo v by p(r, t) and observing that v(r, t) is curl-free. 
Because of the latter we have 

d d 1 o 
— v = —v + -vv 

dt dt 2 



dt 



(P(*)) 



(85) 



which is Ehrenfest's first theorem. 

In case of a force- free particle for which (F) = 0, Eq. ([85]) 
yields 

(p(t)) — const. 

which demonstrates that a free particle exposed to 
Brownian/anti-Brownian stochastic forces does not 
change its momentum on the average, as opposed to a 
particle that moves in a classical "Brownian" environ- 
ment. 

We now want to show that the expectation value of "New- 
ton's modified second law" that we have derived in the 
form of Eq. (f2"5]) . attains, in fact, exactly the form of 
Eq. (j85|) . To this end it is convenient to recast Eq. (l24]) 
as 



Vqp 



rn 



4 too 
u 2 (r,t) 



1 (Yf. 

2 ^ p 



V 2 p 



V ■ u(r,t) 



2 2to 

where we have used Eq.© defining u(r, t). Hence 

j3 



J p{r,t)F QM (r,t) d 3 



(86) 



to 



-p(r,t)Vu 2 (r,t) 



2mo 



p(r,t)Au(r,t) 



d 3 r 



which can be recast as 

d d 
dt dt 



dp dp 
m v — - too v — 
ot ot 



m „ 2 

-— pVv 



where we have added zero in the form of the bracketed 
expression. The real-space integral over this equation 
may be written after reordering 



= 7Tt f v f>(r,t) m v(r,t) d 3 r 



TO 



2v^-- pVv 2 
dt y 



(87) 



The integral on the right-hand side vanishes because of 
the equation of continuity 

^+V-(p«)=0. (88) 

This follows from multiplying this equation by v and per- 
forming a real-space integration. We then have 



v — d 6 r 

v dt 



v u V-{pv)d 3 r (89) 



yZ e " V ' {v v pv)d 3 r + V"e„ / pv-Vv„d 



- J. x> u p v-d 2 r 
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with e v denoting unit vectors. The surface integral as 
been obtained by invoking Gauss' theorem. It vanishes 
since we may assume p\v\ to vanish sufficiently toward 
infinity. Again exploiting the property of v being curl- 
free the second integral on the right-hand side can be 
written 



/ p(v-W)vd 3 r = i / pWv 2 d 3 
Jv 2 J v 



r . 



It follows then from Eq. (f8"9"| that the integral on the right- 
hand side of Eq. ([8T|) is, in fact, equal to zero. Thus we 
have shown that the expectation value of the "vacuum 
force" FQM(r,t) vanishes 



J p{r,t)F QM (r,t) d 3 r = 



which plays also a central role in information theory 
(s. e. g. Garbaczewski^). 



XIV. THE TIME-DEPENDENT SCHRODINGER 
EQUATION IN THE PRESENCE OF AN 
ELECTROMAGNETIC FIELD 

In going through the various steps that led from 
Eq. (f2"3")l ("Newton's modified second law") to the time- 
dependent Schrodinger equation ([70]) one recognizes that 
we implied nowhere that F has to be time- independent. 
Hence one is justified in allowing F in Eq. (|2"3")) to be time- 
dependent and attain the particular form 

F(r, t) = - W cons (r) + e E(r, t) 

+ev(r,t) x B(r,t) (90) 

if the particle under study possesses the charge e and is 
acted upon by an electric field E(r,t) and a magnetic 
field B(r,t). The quantity V cons (r) denotes the poten- 
tial of an additional conservative field (e. g. the gravi- 
tational field) which we include to ensure full generality, 
and v(r,t) is the ensemble average defined by Eq.Q. 
From B = V x A and Faraday's law of induction we 
have V x (E + A) — 0, and hence E + A may be ex- 
pressed as a gradient of a scalar function which we denote 
by -±V el (r,t). Thus 



e E(r, t) = -e A(r, t) - VV ei (r, i) . 



If the magnetic field is switched on, it induces a voltage 
Vr along any circular path C 

V R = £ E md . (r', t)-dr' = ~J B(r', t) ■ d 2 r' 

where C is the rim of the surface A. On multiplying 
this equation by e and observing that e Ei n d. represents 



an additional force that changes the momentum of the 
particle, we obtain 



p(r',t')-dr' = i eE ind .(r',t')-dr' = 
< Jc 

d 



} - I eA(r',t') ■ dr' . 
t Jc 



Integrating this equation from to to t and assuming 
A (r' , to) = we obtain 

- f p(r', t ) - dr' + f p(r', t) ■ dr' = 
Jc Jc 

- <j> eA(r',t) - dr' . 



where 



p(r', t ) • dr' = mo (p v(r' , t ) • dr' = , 
c Jc 



which follows from Eq. (|H]) . Thus 



v(r',t) + — A(r',t) 
m 



■ dr' = V t 



which means that the curl of the integrand vanishes: 



V x 



v(r,t) + A(r,t) 

mo 



(92) 



Consequently, it can be expressed as a gradient of a scalar 
function which we denote by (h/mo) (fi(r, t). Hence we 
arrive at 

e h 
v(r,t) + — A(r,t) = — Vip{r,t). (93) 
mo mo 

which now stands in place of Eq. P5|) . 

We note here only in passing that we have because of 

^(r) = |V>(r)|e^M 

— Up*Vip - tpVip*} = \i>{r)\ 2 V<p . 
2i 

Using Eq.(|93|) one can recast this as 

h e 

[ip*Vip - 4>V%})*] = pv + |V>| 2 A 

2mo Too 



(91) or equivalently 



pv 



2toq 



[ip*P<ip + c.c] 



where P is short-hand for p — e A. After real-space in- 
tegration and an integration by parts one arrives at 



(v) = — I f(r,i)P^(r,t). 
to 



(94) 
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Because of Eq. (l92l) the expression (v-W)v which appears 
in 



The left-hand side of Eq. (|98j) is obviously the time- 
derivative hereof. It will be useful to notice that 



m -j:v(r,t) = m 
at 



— v + (v ■ V)v 

Ot 



F(r,t) + F QP (r,t) (95) 

cannot be replaced by ^Vv 2 any more. Because of the 
generally valid relation 



(a • V) a = V — — ax (V x a) 



and because of Eq. (|92|) we now have 

(v ■ V) v = ^Vt> 2 - v x (V x v) = 

1 e 

-X7v 2 H v x (V x A) 

2 m Q y ' 

Using V x A = B we may recast this as 



(v-V)v 



1. 



-v x B . 



m 



Inserting this result together with Eq. ([9T)]) and -Fqp = 
— VVqp into Eq. ([95|) we notice that the Lorentz- force 
ev(r,t) x B(r,t) drops out in favor of A(r,t), and we 
get 



d_ 

at 



e 

m 



1 



1. 



(v + — A) = — — W - -Vv z 

mo 2 

1 r, 2 ^ A 

+ - Vix 2 Am 

2 2m 



where we have introduced 

V{r,t) = V cona (r) + V e i (r,t). 



(96) 



(97) 



We now multiply Eq. (|96|) by the imaginary unit z and 
subtract Eq. (|60|) which gives in complete analogy to 
Eq.® 



|r — V(ln^/V^) = V(^-^)- (101) 
ot m m ip ot 



We may also use Eq. (|100|) to recast the first expression 
on the right-hand side of Eq. 



V( — u + i v)' 



2 to 2 , 



VlVln^/V^)] 2 



+ A _^ V [A . VlmVA/po")] - ^V(— A) 2 
to to 2 to 

If one observes that 



v[v WWVpd] = -x ^ - I v MW\/po~)] 2 , (102) 

the second expression on the right-hand side of Eq. 
can be written 



i-Ji— V [V • (-u + iv)] 
2 to 



^A_ V [Vln(^/^)] 2 



1 / 1 ft A , 



2 toq toq 



■0 V ■ A ) 

= (V--A)j/>-(A-V)V> 



Hence we obtain 



V(-tt- 



where 



i h 
2 toq 



V[V(-w + ««)] 



■iV 



too "0 



A^ + i V • A -M 2 A • V 

2 too 2 mo 2 too 



2 toq 



(eA) 2 



_5 

at 



u + i(v H — — A) 

too 



- V(-it + i t>) 2 



i 7i 
2 toq 



V [V ■ (-u + iv)] 



W. 



to 



(98) 



We mention here only in passing that Eq. (|60p is equiva- 
lent to Fick's law and is hence not affected by the pres- 
ence of an electromagnetic field as long as Einstein's law 
(|22p remains unchanged which is obvious from his deriva- 
tion. (S. also Fritsche and Haugki^.) 
As in the case without electromagnetic field we ab- 
sorb the two independent scalar informations p(r,t) and 
ip(r,t) into one complex-valued function 



^{r,t)^±^pj^e lip{r ^ 



(99) 



As v is no longer equal to we have now in place 



V(lnV>A/A)) = ~u + i(v+ A). (100) 

m 



of Eq.(pj 
_h 

to 



The right-hand side of his equation may be compactified 
by using the momentum operator (|66p as a convenient 
short-hand notation 



jV(-„ 



IV] 



2 toq 



V [V(-u + iv)] = 



-iV 



1 (p-eA) 2 
mo tp 2 toq 



Inserting this result into Eq. (|98|) which derives from 
Eq.([95| ("Newton's modified second law") and Eq.© 
(= Fick's law) and exploiting the Eqs. fTUTj) and (fTT7^|) 
we arrive at 



iK d -^=H(r,t)i,(r,t) 



(103) 



where 



H(r,t) 



2 toq 



V(r,t) and P = p-eA(r,t) 
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XV. A MODEL FOR NON-MARKOVIAN 
DIFFUSION ILLUSTRATING THE ORIGIN OF 
NON-LOCALITY 

It is instructive to consider a model illustrating 
"conservative diffusion" . The latter is a consequence of 
forming the arithmetic mean of Eqs. fTT)) and (f2"0")) which 
leads to Eq. (|2Tj) . If one were to follow the motion of 
an individual particle, just one member out of the total 
ensemble, one would directly see the effect of stochastic 
forces changing back and forth from "Brownian" to 
"anti-Brownian" with the latter causing a motion 
enhancement after the former have effected a slow down 
of the particle motion. Figure [2] shows three situation of 
the (free) particle which moves within a two-dimensional 
frame where a two-slit diaphragm has been inserted on 
the left-hand side. The "walls" of the frame are assumed 




FIG. 2: Trajectory of the test particle undergoing reversible 
scatterings 

elastically reflecting. The stochastic forces acting on 
the particle are simulated by a two-dimensional gas 
of N identical point masses (N ^> 1) that interact 
via Lennard-Jones pair-potentials with each other and 
with the particle under study as well. The latter will 
henceforth be referred to as "test particle". It is this 
situation which the original derivation of Einstein's 
law (|22p refers to where the motion of the test particle 
is described by a Langevin equation into which the 
embedding of the particle enters through a stochastic 
force. (The practical calculations have been performed 
with slightly modified Lennard-Jones potentials that 
were truncated at twice the average particle distance.) 
In our model the particle motion of the embedding gas 
results from a molecular dynamics simulation which 
one starts by first keeping the test particle fixed at the 
point ta and letting the N gas particles start from some 
corner of the frame with equal absolute values of their 
momenta. Thereby one defines a certain value of their 



total kinetic energy E^f. After a short simulation time 
the gas particles are uniformly distributed within the 
frame and their distribution in the momentum space 
has become Maxwellian. The latter is associated with 
a certain temperature such that the thermodynamical 
expectation value of the kinetic energy equals El?*!. It 
is this temperature which finally shows up in Einstein's 
law (|22p. After the embedding gas has "thermalized" 
one imparts a certain momentum p on the test particle 
and continues the molecular dynamics simulation with 
the test particle now included. As indicated in the upper 
panel of the figure, it performs an irregular (Brownian) 
motion and loses momentum to the embedding gas 
whose particles are not shown in the figure. We have 
chosen the starting point ta such that the particle moves 
through the upper slit of the diaphragm and reaches the 
point tb after a simulation time At of the order of r 
which is the time constant of a freely moving particle 
in a gaseous medium with friction. We now look for a 
point rc further to the right in the forward direction of 
the test particle (s. panel in the middle of the figure). 
At this point we impart a momentum — p on the particle 
(after thermalization of the embedding gas), i.e. just the 
reverse of the momentum at ta ■ 

The point rc is chosen such that the trajectory ends - 
again after an identical simulation time of At seconds 
- at point tb- At this point the test particle has lost 
its original momentum ~p almost completely. If one 
now turns the velocities of all particles around by 180° 
and starts the simulation again with the time running 
forward as before, the test particle continues its motion 
from point tb and moves exactly along the trajectory it 
had formerly followed in the opposite direction coming 
from rc- When it has reached rc again, it has regained 
the previously lost momentum, but this time with 
the sign reversed. Hence, in moving from r a to rc 
the particle undergoes scattering processes that are 
in alternating succession Brownian and anti-Brownian 
within a time interval of the order r. Thereby the 
average momentum of the particle is conserved. This 
is illustrated in the third panel (bottom). A striking 
feature of the momentum reconstruction by the above 
scattering processes is the occurrence of non-locality. 
This can be demonstrated by repeating the procedure 
that led to the trajectory portion from rc to rs with 
a crucial modification: If one closes the lower slit of 
the diaphragm and starts then with the same posi- 
tion/velocity configuration of all particles as before, 
the trajectory of the test particle evolves now differently 
and does no longer join the previously generated trajec- 
tory portion at rs- This is what the molecular dynamics 
simulation clearly yields. On the other hand, this is to 
be expected anyway because every momentum transfer 
from the test particle to the gas spreads with sound 
velocity throughout the entire structure and probes 
the change that has been introduced. The stochastic 
forces acting on the test particle are modified by such a 
change when these sound waves are reflected back on the 
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particle. If one wants the modified trajectory to join the 
first trajectory portion at rg again, one has to choose a 
different starting point r' c . Once the test particle has 
arrived at r^, one inverts all the velocities as before, 
and the particle will now recover the momentum p on 
its modified trajectory toward r' c . Note: this change 
in the course of the particle motion results just from 
closing the lower slit although the particle definitely 
traverses the upper slit. One is tempted to surmise 
that this mechanism of probing the environment "in real 
life" as the particle exchanges temporarily momentum 
with the vacuum, occurs at light velocity. The latter 
would impose a limit on the distance beyond which 
a previously passed potential structure can no longer 
affect the evolution of the particle's trajectory at its 
current position. 

If one were dealing with Brownian scattering only, 
the succession of scattering events could be classed as 
"Markovian" . (Shorthand definition: given the presence, 
future and past are independent.) However, the overall 
character of the combined Brownian/anti-Brownian 
scattering processes is obviously non-Markovian. It is 
true that the particle has almost completely lost its 
memory of its original momentum when it arrives at 
rs, but its future time evolution while moving toward 
rc reconstructs, so to speak, past scattering events. 
The particle's momentum p(0) when it is at r a 1 and its 
momentum p(t) at rc are strongly correlated. 
This does not apply to the positions ta(0) and rc(t): if 
one repeats the experiment and lets the particle start at 
ta with the same momentum p(0) as before, but with 
the thermalization process of the embedding gas started 
some time interval earlier, the particle's trajectory 
will now be different and lead to a point different 
from rc though it regains its original momentum after 
(approximately) the same traveling time. 
Obviously, the non-Markovian (reversible) character of 
particle motion which results from such a combination of 
scattering processes can only show up on a coarse grain 
time scale which is the crucial assumption underlying our 
derivation of Newton's modified second laws Eqs. (f2"5|) 
and §5$ together with flU}. 



XVI. OPERATORS AND COMMUTATORS 

An important advantage of our approach may be seen 
in the derivability of Hermitian operators which in stan- 
dard quantum mechanics can merely be obtained from 
educated guessing employing Jordan's replacement rules. 
In Section [X] we have already derived the momentum op- 
erator 

p = —i HV 

exploiting our expression (|4"S"|) for v(r,t) and j(r,t) = 
p{r,t)v{r,t). The same arguments used in deriving p 



apply to the angular momentum operator L which oc- 
curs on forming the expectation value of the angular mo- 
mentum of a particle with respect to a center located at 
r = 0. This expectation value (L) is primarily defined as 
a real-space integral over the angular momentum density 
r x moj: 

(L(t)) = J r x m j(r,t)d 3 r (104) 

If one here inserts j from Eq. lfTTj) . integrates by parts 
and requires V , ( r j t) to vanish sufficiently toward infinity, 
the result may be written 

(L(t)) = J \p*{r,t)(r x p)i;(r,t)d 3 r (105) 

which justifies terming L = r x p "angular momentum 
operator" . 

The kinetic energy of an individual particle, labeled by 
the index j, is defined as the work performed on that 
particle by the external force F in accelerating it from 
zero velocity at time t = to its velocity v C j(t) at time 
i, which yields 

Forming the ensemble average according to Eq.Q one 
obtains E kin = ^v^(r,t). Thus, the density of the 
kinetic energy is given by 

efe in (r, t) = ^Y p(r, t) [v c (r, t)f . 

In the two subsystems "B" and "A" we are considering 
the velocity v is always the same, whereas the convective 
velocity v c is different, and therefore we distinguish 
from v£ and form the ensemble average over the two 
subensembles: 

e kin (r , t) = m \ ( [vf (r, t)] 2 + [ttf (r , *)] 2 )(106) 

According to Eq.((6|) which still refers to the "B" -system, 
we have 

— v — u and therefore = v + u . 
Consequently Eq. (|106p may be cast 

e hin .{T,t)=m ^- {[v(r,t)] 2 + [u(r,t)} 2 ) .(107) 
From Eq. ([62|) we have 

fi 1 r, , 
— u + i v = Vip . 

mo tp 

The modulus square of this equation times m p/2 is 
equal to the right-hand side of Eq. (|107[) . that is 

h 2 

efcm.M) = - |VV(r,i)| 2 . 

2 ran 
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Taking the real-space integral of this expression one ob- 
tains the kinetic energy 



E km = (T(t)) = 



2 mo 



Vip*(r,t) ■ Vip(r,t) d 3 r (108) 



which by employing Green's theorem may be given the 
familiar form 



v 2 mo 



Vif)*(r,t) ■ V^(r,t) d 3 r = 



2 toq 



ip(r,t)d 3 r , 



and hence 



(T(t)>= / ^(r,i)-^-V(r,i)d 3 r, 
Jv 2m o 

which justifies terming p 2 /2mo "kinetic energy opera- 
tor". 

In practical calculations one often benefits from the fact 
that Ekin may alternatively be cast as in Eq. p08[) where 
the integrand is real-valued and may immediately be in- 
terpreted as "kinetic energy density" . 

The statistical operator is a particular example of deriv- 
ability from a simple concept. We confine ourselves here 
to the case of a quantum mechanical system of a bound 
particle in contact with a heat bath of temperature T. In 
a stationary state the latter constantly exchanges energy 
with the system, in the simplest case photons. Hence, 
the wave function of that system cannot be one of its 
eigenstates any more, but rather represents a solution 
to the time-dependent Schrodinger equation and can be 
expanded in terms of eigenfunctions VVi( r ) 



VXM) = J! c n (t)ip n (r)e 



(109) 



where E n denotes eigenvalues of the unperturbed one- 
particle Hamiltonian H . 

To make the external system classifiable as a heat bath, 
the time-averaged coupling energy of the two systems 
must be negligibly small compared to the difference 
E n i — E n of any two eigenvalues. The particle's ther- 
modynamical expectation value of its energy (indicated 
by double brackets) is given by 



I + t r 



((H)) = U = 
ij*(r,t')Hi;(r,t')d 3 r 



dt' 



(110) 



where r (not to be confused with the slow-down time in 
Section HV|) has to be chosen sufficiently large such that 
U does not depend on t any more. Quantities that derive 
from U like the specific heat, are only defined as time- 
averages of this kind. 
Inserting Eq. (|109p into (|110|) we obtain 



U 



\c n (t')\ 2 dt'} 



(111) 



where the expression in curly brackets may be interpreted 
as the relative frequency of the system of being in the n- 
th eigenstate. 

Straight-forward thermodynamics yields for a system 
that possesses energy levels E n 



U^Ej-e-^, f,= * 



k B T 



(112) 



where 



Thus, we have from Eq. ljlll 



' ft+T \ Cn (t')\ 2 dt' = -e^ E - 



If one defines a statistical operator 



_ 1 
p=-e 
a 



-PH 



Eq. (|112j) can alternatively be cast as 

E/ = X>«IP%n) =Tr(pH). 

n 

Commutation rules for the operators apply when 
the potential V(r) in the time-independent Schrodinger 
equation ([50)) possesses a certain symmetry. If V(r) is 
spherically symmetric, for example, one verifies simply 
by performing partial differentiations that 

HL 2 - £ 2 i?) Vn(r-) = [H, L 2 } Vv,(r) = 

and similarly 

[H,L z ]i> n {r) = 
if VVi( r ) is an eigenfunction of H . 

If one is dealing with some operator A which represents 
just some analytical expression in r and p, the time de- 
pendence of its expectation value (A) can be determined 
by employing the time-dependent Schrodinger equation 
which gives 



dt 



ip*(r,t)Aip(r,t) d 3 r = 



r(r,t)-{H,A]^r,t)d 3 r., 



in short-hand notation 

Commutation rules of the above kind, again in short- 
hand notation 

[H,L 2 ]=Q; [£,£,] =0, 
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similarly 

[H,p\ = if V(r) = const., 

but also 

\pj,x k ] = -5j k where j = 1, 2, 3 ; k = 1, 2, 3 

constitute fundamental elements of standard quantum 
mechanics and are discussed as pivotal in the context 
of measurement. From our point of view they are just 
byproducts of the Schrodinger equation and do not con- 
tain any more physics than has already gone into the 
derivation of the Schrodinger equation. In practice it 
is impossible to find quantum systems where eigenval- 
ues of H and L z , for example, can be measured simul- 
taneously although there is a widespread belief to the 
contrary. It is not even possible, for example, to mea- 
sure the eigenvalues of H for a hydrogen atom which - in 
clamped proton approximation - represents the archety- 
pal one-particle system and the starting point of quantum 
mechanics. The lines one observes in its discrete optical 
spectrum refer to eigenvalue differences and possess - dif- 
ferent from true eigenvalues - a natural line width which 
goes to zero only in the hypothetical case of zero radia- 
tion coupling, that is when the lines cannot be observed 
any more. 

There is a remark by Wigner 40 which reveals exactly that 
lack of stringency and consistency in the foundation of or- 
thodox quantum mechanics: " All these are concrete and 
clearly demonstrated limitations on the measurability of 
operators. They should not obscure the other, perhaps 
even more fundamental weakness of the standard theory, 
that it postulates the measurability of operators but does 
not give directions as to how the measurement should be 
carried out. " 



XVII. COLLAPS OF THE WAVE FUNCTION 
AND THE NODE PROBLEM 

A vital point of the Copenhagen interpretation consists 
in the notion that the wave function of a stationary one- 
particle state collapses on performing a measurement on 
the position of the particle, for example. Within our ap- 
proach a phenomenon of this kind cannot occur. First of 
all, in our view "measurement" is not a process of some- 
thing foreign intruding the realm of quantum mechanics 
but is rather a part of it. If one calculates, for example, 
the time-independent wave function ip(r) for a stationary 
situation where electrons in a diffraction chamber leave a 
tunneling cathode, sufficiently far behind each other, run 
through a two-slit diaphragm and finally hit a fluores- 
cent screen, ^(r)! 2 will display the familiar diffraction 
pattern behind the diaphragm and in particular on the 
screen. But clearly, the structure of this pattern reflects 
the distribution of the entire ensemble of electrons that 
leave the cathode, and a particular electron, that hits the 



screen some place, is only one member out of this ensem- 
ble. Hence its capture on the screen does not destroy the 
properties of the ensemble. The electron capture by an 
atom of the screen constitutes a process that has only 
marginally to do with the diffraction state in that the 
latter determines the probability of the electron being 
at that particular atom. Otherwise the capture process 
is governed by the time dependent Schrodinger equation 
and the perturbation caused by the electromagnetic field 
of the outgoing photon. All this is completely indepen- 
dent of the possible presence of an "observer" who might 
see that photon. 

Despite deceptive similarities the situation becomes con- 
ceptually different when one replaces the tunneling tip in 
the otherwise unchanged diffraction chamber by a light 
source that emits, again in sufficiently large time inter- 
vals, photons of the same wave length as the previously 
considered electrons. Since the space-time structure of 
the wave (in principle oc cos[fc ■ r — cot]) with which each 
photon is associated, is not defined as the property of an 
ensemble of mechanical objects but rather by classical 
electrodynamics (Maxwell's equations), it will, in fact, 
disappear on the disappearance of the photon in ques- 
tion. There are a couple of properties by which photons 
differ crucially from massive particles: They move always 
at light velocity along straight lines in vacuo and the 
associated waves are vector-valued functions. By con- 
trast, the de Broglie waves of massive particles are in 
general complex- valued functions, and the average ve- 
locity of the particles is given by the gradient of the 
functions' phase. Their interaction with other particles 
and parts of an experimental setup is described by the 
Schrodinger equation and the potentials therein. On the 
other hand, the interaction of photons with polarizers, 
mirrors, quaterwave plates, filters etc. is governed by 
classical electrodynamics. Malus' law, for example, con- 
stitutes a law of classical optics. Because of these rather 
fundamental differences any analysis of photon correla- 
tion experiments, for example, should critically be scru- 
tinized whether a transfer to analogous experiments with 
massive particles is truly justified. In Section IXX VI 1 1 1 we 
shall draw on the familiar example of the Stcrn-Gcrlach 
experiment to demonstrate that the selection mechanism 
for up-spin and down-spin particles in the Stern-Gcrlach 
magnet has nothing to do with the mechanism separating 
horizontally and vertically polarized photons in a polar- 
izing beam spitter. 

According to Miclnik and Tengstrand 3 - excited station- 
ary states appear to pose a serious problem in that ip{r) 
possesses nodal surfaces at which the normal derivative 
-g^pix) vanishes but the normal component of the os- 
motic velocity 



becomes formally infinite. Moreover, at surfaces across 
which p(r) attains a maximum, ^p(f) vanishes as well, 
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but u n (r) becomes now zero. If ip(r) is real-valued then 
v(r) vanishes everywhere, and therefore we have on such 
surfaces with maximum probability density 

v = u n = . 

In the 2s-state of a hydrogen electron, for example, one 
has a spherical surface of this kind. Hence, it seems that 
this sphere separates two regions of space that are mutu- 
ally inaccessible for the electron. But the above velocities 
are only ensemble averages or - in the spirit of the defini- 
tion ([5]) - averages of non- vanishing velocities v(ti), u n (fi) 
of different directions over a sufficiently long time T. 

As for |it„(r)| going to infinity as one crosses a nodal sur- 
face of tp(r), one has to keep in mind that stationary ex- 
cited states (excited eigenstates) are highly fictional and 
do actually not exist in nature. Because of AE At as fi 
and AE = for an eigenstate it would take an infi- 
nite time to prepare them. Hence, truly existing excited 
states do not possess nodal surfaces where ip(r) vanishes 
exactly. But even if one would allow them to exist, the 
kinetic energy density ^p(r)u 2 (r) — V?/>( r )| 2 re- 
mains finite and hence ensures a physically meaningful 
behavior even for this idealized situation. 



XVIII. THE FEYNMAN PATH INTEGRAL 

As our concept builds on the existence of particle tra- 
jectories one might surmise that there should be some 
affinity to Feynman's path integral method^ which also 
relates to possible paths a particle might take. We shall 
outline that there is neither any formal kinship nor does 
Feynman name any cause for the possible occurrence of 
non-classical trajectories. In so doing we limit ourselves, 
as Feynman in his article, to the one-dimensional case 
of a particle that moves non-relativistically in a poten- 
tial V(x). Feynman's considerations are based on two 
hypotheses that may be summarized by stating that the 
wave function tp(x,t + At) of the particle at some point 
x and time t + At is connected with the wave function 
ip(x — a,t) at a previous point x — a and earlier time 
t by an integral equation similar to the Smoluchowski 
equation (|215[) of the ensuing section, viz. 



V(x) along a trajectory from a point x — a to x within 
an infinitesimally small time span At. Hence 



S(x, x — a, t, At) = 

nt+At 



Min. 



^a*-V(x-a(t>)) 



dt' 



where 



L(a(t),a(t)) = ^a 2 -V(x-a(t)) 
denotes the Lagrangean. 

As At is infinitesimally small S(x,x — a,t,At) may be 
approximated 



S = At 
Hence one has 



m / a y 
2 \At) 



V(x) 



F = 



' '"0 -2 i V(x) At 

g 2h At u . g h 



The first exponential oscillates rapidly as a function of a 
because of the prefactor I /At in the exponent, whereas, 
by comparison, ?p(x—a, t) may be assumed slowly varying 
as a function of a. The value of the integral in Eg. (|1 13[) 
depends therefore only on a small interval of a around the 
point x. Within this interval ijj(x—a, t) may be expanded 
as 

/ s 1 1 \ dil) 1 d 2 ip o 

ip(x -a,t) = ip(x, t) - — a + - — a . 

If one inserts this into Eq. (|113[) , observes 



f 

A 



e 5 ™ da = 1 : 



(note that this equation defines A\), further 



e 2K ff2 a da = , 



ip(x, t + At) = J ip(x - a, t) F(x, x-a,t, At) d<r(113) 

where F(x,x — a,t,At) is the function that brings in 
classical mechanics. It is defined as 

1 i 



F(x,x - a,t,At) = — e %S{x,x-<r,t,At) 



where 



and uses 



e^TAt a 2 da = At 



A 



j vjx) At iV(x)At 

R Ri 1 



one obtains 



2nhiAt \ 2 
mo 



Here S(x, x — a,t, At) denotes Hamilton's first principle 
function for a particle moving classically in a potential 



ip(x, t + At)= ip(x, t)\^l-- V(x) Atj + 

1 d?ijj ih / iV(x)Af 

2 dx 2 mo \ Ti 
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Multiplying this equation by j& and letting At tend to 
zero one arrives at the time dependent Schrodinger equa- 
tion 

h 2 d 2 

2 mo dx 2 



d 



V(x) 



■>jj(x, t) 



Though the Schrodinger equation is obviously recovered 
following this line of argument, it remains unclear why 
Hamilton's classical first principle function should appear 
in the exponent of F(x, x — a, t, At). Feynman's consider- 
ations lean closely on arguments of measurement typical 
of the Copenhagen school of thought, cast into an ax- 
iomatic framework notably by v. Neumann--. But tj}(x, t) 
may, for example, describe the motion of a harmonic os- 
cillator in the absence of any measurement. In fact, if 
one were to perform a measurement on the harmonic os- 
cillator the Schrodinger equation would contain a per- 
turbative extra term that would give rise to a different 
wave function. Clearly, as already stated in Sectiorfl] the 
probabilistic character of ip(x,t) does not originate from 
indeterminacies caused by the process of measurement. 
The complex- valuedness of the wave function in the form 
ip(x,i) = \ip(x,t)\ e 1 ^^'^ comes about by incorporating 
two autonomous real-valued informations: the probabil- 
ity density \ip(x,t)\ 2 of the particle being at x and time t 
and the ensemble average v(x, t) = <-p(x, t) of its ve- 

locity. For that reason our derivation of the Schrodinger 
equation requires two Smoluchowski equations for the 
real- valued functions p(r,t) and v(r,t) instead of Feyn- 
man's single Eq. (|113p . We believe, therefore, that in our 
derivation the connection to classical mechanics becomes 
definitely more transparent and convincing. 



XIX. THE TIME-DEPENDENT N- PARTICLE 
SCHRODINGER EQUATION 

So far we have merely been concerned with a single 
particle whose stochastic behavior was described by re- 
garding it as a member of N identically prepared, but 
statistically independent one-particle systems under the 
supposition that N be sufficiently large. To avoid con- 
fusion we shall henceforth rename that number by Af. 
Instead of a single particle we now consider N particles 
that interact via pair- forces. Each of these particles is in- 
dividually a member of Af statistically independent one- 
particle systems where the N — 1 remaining particles ap- 
pear at fixed positions r 2 ,r 2 , . . .r^ if the particle under 
consideration, picked at will, just happens to be "number 
1". The considerations of Sections IXl and IXIVI carry over 
to this iV-particle system. To see that one simply has 
to replace the 3-dimensional real-space of the single par- 
ticle discussed as yet by a 3-ZV-dimensional space where 
the N particles appear as one point again. Instead of the 
probability density p(r,t) one is now dealing with 

N 



p(n,r 2 ,...r N ,t) = p(r ,t) 



where / p(r N , t) d 6 r x d J r 2 . . . <?r N = 1 (114) 



and 

N 3 

r N = (r 1 ,r 2 , ■ ..r N ) = ^ Xj k ej k 

j=i k=i 

with j — 1, 2, ... N numbering the particles and Xj k de- 
noting Cartesian coordinates which are associated with 
orthogonal unit vectors ej f. . The quantities , u N and 
v N are defined analogously. 
Instead of <p(r,t) we now have 4>{r N ,t), Thus 



m 



(115) 



Correspondingly, the 3iV-dimensional osmotic velocity 
has the form 



u N {r l \t) 



2 mo 



V JV ln[p(rV)/Po], (H6) 



and hence we have similar to the single-particle case 



du 



N 



7 N 



dt 



d_ 

2 mo " dt 
h _ 



[lnp/po] = 



iV 



1 dp 
~p dt 



(117) 
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Invoking the equation of continuity 

|+ V^Q =0 

and using the definition (|116p . Eq. (|117|) can be cast as 
du N h 



Ot 



2 mo 



V w [(V" -v N )-{u N •/)]. (118) 



In the following we first confine ourselves to time- 
independent conservative forces which - in the spirit of 
our notation - may be written 



?N 



N 3 

j=i fe=i 



where 



■(r j ) = - 



dx 



3 



V ext .(rj) 



with Vext. ( r ) denoting an external potential. Hence F t 
may alternatively be written 

F» t Xr N ) = -V N V ext .(ri,r 2 ,...r N ) 



N 



where 



N 



V ext .(ri,r 2 , ■ ■ -r N ) = y^ y V e xt.( r j) 
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The force exerted on the j-th particle due to pair- 
interaction with the N — 1 remaining particles is given 
by 



Here Hq denotes the "free Hamiltonian" 



F inter {rj) = 



d 



N 



dx 



ok 



Hq = Hj where Hj 

3=1 



2 nio 



+ Vext.{rj) 



.(122) 



where V(\rj — Vi\) denotes the interaction potential. The 
generalized total force in the 37V-dimensional space may 
therefore be cast as 

F N (r N ) = -V N V(r u r 2 ,...r N ), 
where V(ri,r 2 , ■ . .rjv) is defined by 



V(r 1 ,r 2 , ■ ..r N ) = ^ V ext . fa) 



Because of Eq. (|115l) the phase of the wave function 
may still depend on time when p and v N are time- 
independent: 

cb(r N ,t) = Mr N ) + f(t). 
Thus we have in this case 

*(rV) = ^o(r N )e- tfW ; 
* (r JV )=± v /p(r^)exp[z ^(r w )] 

which on insertion into Eq. (112ip yields 



N N 



I^^y(| r ^_ ri |).(119) Hy a (r N )=hf*o(r N ) %f = const. = E 



3=1 

Newton's modified second law (9.2,1) hence attains the 
form 



E 



dv 



N 



at 



— V+\{v N f- l -{u N f- 
mo 2 2 



■ V N ■ u N 



whereby Eq. (|121[) becomes the time-independent 
Schrodinger equation 



H^ Q (r N ) = E^ (r N ) 



(123) 



2 mo 



(120) 



As in the one-particle case the two scalar functions 
p N (r N ,t) and 4>(r N ,t) can be absorbed into a complex- 
valued function ^fa, r%, . . .tn, t) defined by 



*(ri,r 2 , ...r N ,t) = ±y / p(r 1 ,r 2 , . ..r N ,t) x 

cxp [i(j)(r 1 ,r 2 , . ..r N ,t)} . 

This is equivalent to 

-u N (r N ,t) + iv N (r N : t) = —V N (hi*(r N ,t)/^pE) 



XX. STATES OF IDENTICAL PARTICLES AND 
ENTANGLEMENT 

If the particles are non-interacting, one would naively 
expect their motions to be completely uncorrelated which 
means 

N 

p( ri ,r 2 ,...r N ,t) = ]]_Pj{r j ,t), (124) 

3=1 



and 



n 



m * 

which is the analogue to Eq. (9.9,1), and we obtain ac- 
cordingly 



(ri,r 2 , ...r N ,t) = ^2<Pjfa,t) . (125) 

3 = 1 



-I, 

dt 



-u N + iv N ) 



7 N 



n i a* 

m * dt 



If we here insert Eqs.dHl) andCOP for 2*1 an d 
and proceed exactly as in the single-particle case we ar- 
rive at the ./V-particle Schrodinger equation 



*(ri,r 2 , . ..r N ,t) 



=H 



ih—y(r u r 2 ,...r N ,t). (121) 



In that case Eq. (|115p attains the form 

(v 1 (r 1 ,t),v 2 (r 2 ,t) 7 . ..v N (r N ,t)) = 

— (Vnpi(ri,t),V 2 tp 2 (r 2 ,t),...V N p N (r N ,t)) . 
m 

Likewise, Eq. (|116p becomes 

(tti(ri,t),« 2 (»"2,t), ■ ■■u N (r N ,t)) = 

-7T— (Viln[/3i(ri,t)//j i], 
2 mo 

V 2 m[ / 9 2 (r"2,£)/po2], ■ . . Vjv[lnpAr(rjv,i)/poJv]) ■ 

Newton's modified second law (|120[) decomposes accord- 
ingly into N analogous equations for single particles, as 
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has to be expected. Each of these equations can be 
subjected to a Madelung transform which yields time- 
dependent one-particle Schrodinger equations solved by 
one-particle wave functions tpj (rj ,t). If one multiplies 

Hj(rj) ipjir^t) = ih ^ ^(r i( t) (126) 



by Y\i=i 4 , i(fi,t) one obtains 



N 



N 



Hji r j) II ^(n,t) = Y[ipi(ri,t)ih — ipj(ri,t) 



i = \ 



which on forming the sum 53j=i yields, in fact, 

d 

H ^(r l7 r 2 , . . .r N ,t) = ih — *(n,r 2 , . ..r N ,t) 

N 

where Y&(n.,r 2 , • • . r N , t) = JJ ipj{r h t) . (127) 

Hence, the above time-dependent TV-particle Schrodinger 
equation is solved by the product of individually time- 
dependent wave functions i/jj (rj , t) . 

Obviously, the density (|124[) that results from this wave 
function is not invariant against interchange of any two 
particles if they are in different states, say ip kn {f k ,t) and 
q Pi m ( r h t ) where k n ^l m . 

It is not exactly physical wisdom but rather firm belief 
that even non-interacting massive particles, though non- 
existing in nature, do not perform an uncorrelated mo- 
tion and can, therefore, not be described by the wave 
function (|127p . This belief is based on the idea that 
the particles cannot be tracked individually as they move 
(contrary to classical particles) because the uncertainty 
relation "forbids" the existence of trajectories. Our ap- 
proach to the many-particle problem is characterized 
by the plausible assumption that each particle can be 
identified any time by an affix if it has been assigned 
to a certain number at some chosen instant since each 
particle follows an individual trajectory. The quantity 
v i{ r ii r 2, ■■■ fN, t), for example, represents the average 
over all particle velocities at r\ and time t of the ensem- 
ble associated with particle "number 1" . In forming this 
average the positions r 2 , r 3 , . . . of the N — 1 remain- 
ing particles are kept fixed, that is, the average results 
from the entire set of "number 1" -trajectories that occur 
in the "number 1" -ensemble while the particles "number 
2, 3 „JV" are at fixed positions. Clearly, if one or more 
of those particles are kept at different positions and if 
all particles interact, V\ will in general be different at r\ 
and time t. Hence in our view there is no extra quantum 
phenomenon of indiscernibility. As in classical mechanics 
it is entirely sufficient to characterize identical particles 
merely by their property of having the same mass and 
charge. 

If one insists, however, on "quantum indiscernibility" also 



for non-interacting particles, that is, on the invariance of 
p( r i, f*2> • • • f"zv) against interchange of any two particles, 
one has to replace (|127l) with a renormalized linear com- 
bination of all Nl products that differ in the interchange 
of two particles 

*(ri,r 2 , . ..r N ) = 

N\ N 

-= J2 (±V P Hk, I) I] VS- fo) • (128) 
viv. p=1 j=1 

where P(k, I) is the permutation operator exchanging the 
particle referring to j = k with that for j = I, and P num- 
bers the permutations. 

If the particles interact and are bound in an ex- 
ternal potential or move in a parallelepiped where 
^(ri,V2, ■ ■ -Tjv) is subjected to periodic boundary con- 
ditions, each particle is constantly scattered, and hence 
the probability of some particle, say "number k" , being 
within an elementary volume A 3 r around r is given by: 

P(r) = J \V(r 1 ,...r k - 1 ,r,r k +i ■ ■ ■ r N )\ 2 d 3 n . . . 

d 3 r k _ 1 d 3 r k+1 ...d 3 r N . 

Indiscernibility means that P(r) is the same for any 
particle one picks, that is, each particle appears at r with 
the same probability. Hence we have 

p(r) = NP{r) 

with p(r) A 3 r denoting the probability of any of the N 
electrons being in A 3 r. 

The function p{r\, r 2 , . . . rjy) is now naturally invariant 
against interchange of any two particles. 
An important property of particles is their spin which will 
be discussed farther below in this article. In the present 
context it may be sufficient to introduce 

x = (r, a) 

as a generalized particle coordinate where a = ±1 de- 
notes its discrete spin coordinate and refers to paral- 
lel or anti-parallel orientation with respect to a global 
axis. The wave function (|128[) for non-interacting parti- 
cles then takes the form 



*(a3i,x 2 , . ..x N ) 



N ! 



N 



^(±i) p p(M)n^w- (129) 



p=i 



The alternative in the sign under the sum is related to 
the two fundamentally different species of particles: The 
plus sign in (±l) p characterizes bosons, the minus sign 
fermions. Hence the latter are associated with a wave 
function that changes sign on interchanging any two par- 
ticles. This property persists when *S>(xi,x 2 , ■ ■ ■ x N ) de- 
scribes N interacting fermions. Antisymmetry of the 
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wave function gives rise to a peculiar behavior of the 
so-called pair-density 

/ Dot 
p 2 {x,x ) = 

N(N-l) [ \y(x,x',x 3 ,...x N )\ 2 d 4 x 3 ...d 4 x N 



the two electrons is given by 



where 



Obviously 



,d 3 r 



. . .x„,x v+1 , . ..x N ) ee ^{x x , . . .x u+ i,x v , . ..x N ) 

if x v = x v+ i . 

On the other hand, ^ is required to change sign on in- 
terchanging two particles, and hence the above equation 
can only hold if equals zero if the coordinates of any 
two particles are equal. Thus 

P2(x,x') = if x'=x. 

This indicates the occurrence of the so-called Fermi-hole 
which is absent in bose-particle systems. 
The form of the wave function (|129|) may be cast 
as a determinant, named after J. C. Slater. In so- 
called EPRB-experiments (EPRB=Einstein, Podolsky, 
Rosen^, Bohm^i) which were originally devised to test 
possible correlations between two macroscopically distant 
fermions in a singlet state, the associated wave function is 
just a 2 x 2 determinant. The respective two one-particle 
states are in this context commonly referred to as "en- 
tangled states" . 

The requirement of antisymmetry, which is equivalent to 
the Pauli exclusion principle, is a strong subsidiary con- 
dition in solving the Schrodinger equation (|123|1 . Wave 
functions associated with fermions constitute only a small 
subset of the set of functions that satisfy the Schrodinger 
equation (|123[) . 

It should clearly be stated that the antisymmetry of 
the wave function is definitely not a consequence of our 
stochastic approach, but rather has to be required as an 
additional property, as in standard quantum mechanics. 

The derivation of the time-dependent Schrodinger equa- 
tion (|12ip can again be extended to the case where the 
particles move in an electromagnetic field. The exter- 
nal potential becomes time-dependent then (V ex t.{r) — > 
V{r,t)) and pj has to replaced with Pj(r,t) = jij — 
eA(rj,t). 



XXI. A BORDERLINE CASE OF 
ENTANGLEMENT 

We consider a hydrogen molecule whose nuclei are lo- 
cated at Ra and Rb, respectively. The Hamiltonian of 



H 



E 

fe=i 



(-ihV k -eA(r k ,t)) 2 
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+ 



E 



V(r k ) 
1 



(130) 



where 



V(r) 



47re |r-fi A | 4tt e |r - R B \ 



and eo denotes the electric constant. 
We first assume that there is no external field (A(r, t) = 
0) and that the 2-electron wave function has for large 
proton- proton separation, that is when Rab = \Ra — 
Rb \ 3>Bohr radius, still the entangled form of a singlet 
state dictated by the Pauli principle 

*(r-i,r 2 ) = 

[±{n , T) ® ±{r 2 , 1) - ±(r 2 , 1) ® ±{ri , ?)] (131) 

where 

±{r,a) = [a a (R AB )(p A (r)+b <T (R AB )<PB(r)]x(<r)(l32) 
and 

ct=T(I); al + bl = l; <p A/B (r) = <p(r - R A/B ) 
with the property 



J \^(r-R A/ B)\ 2 d 3 r = l 



Here the integrand denotes the electronic ls-orbital of a 
single hydrogen atom, and PA/ B ( r ,t) = \ f( r — Ra/b)\ 2 
is the associated probability density. Furthermore, the 
unit spinors x(c) have the property 

xV)x(<x)=<W 

Under the supposition that Rab = \Ra — Rb\ is suffi- 
ciently large, say 10 cm or even larger, the expectation 
value (H) attains a minimum for either 

(a T - 1 ,a x -> 0) — (6 T -» 0,64 -» 1) "case 1" 



or 



(«l 



, a 



1) 



(&t — 1 , 6 A 



0) 



case r 



For both cases (H) yields the correct value, viz. -2 Ryd, 
as has to be expected for two hydrogen atoms, each of 
which possesses the energy -1 Ryd. Since the spin struc- 
ture does not reflect the symmetry of the potential, one 
forms a symmetry-adapted linear combination 



*i(ri,r 2 ) 



V2 



[ty(ri,r 2 ) +* r (n,r 2 )] , 
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where 



and 



VA(ri)x(T) ^A(r 2 )x(T) 
vbOOxOL) <fB(r 2 )x(l) 



Pu(n)x(T) ¥>B(r 2 )x(T) 
^(n)x(|) ^(r 2 )x(|) 



The two 2-electron functions are associated with the same 
energy which hence applies to ^i(ri,r 2 ) as well. As a 
consequence of the symmetry of 5 , i(r 1 ,r 2 ) in r% and r 2 
we have 



P(ri) = 

and 

p{r 2 ) 
Moreover 



|*(ri,r2)rcfra = p A (ri) + Ps(r x ) 



|*(ri,r 2 )rti d ri = p^(r 2 ) + p B (r 2 ) ■ 



PA/B(ri/ 2 ) d 3 r 1/2 = - and 



PA/Bir^d 6 ^ + / p A /B(r 2 ) d A r 2 = 1 . (133) 



That means that each electron appears in each of the 
atoms (A and B) with the same probability. This has 
rather implausible consequences if one exposes, for ex- 
ample, one of the atoms (say A) to a Laser puis of fre- 
quency uj. Now A(r, t) is no longer zero. The associate 
perturbation operator has the form 

V pert urb{ri,r 2 ,t) = 

Y?k=i A ( r k,t) ■ V fc if ri,r 2 in or near atom A 
else 

which promotes the 2-electron system to an excited state 
*/(ri,r 2 ) = [^ f) ( ri ,r 2 ) + ^P(r u r 2 )] 



where 



* z (/) (ri,r 2 ) 



V2 



and 



^/\r lt r 2 ) = -±= 



V ( l\n)x(V ^ /) (r 2 )x(T) 
yB(ri)xU) ¥>B(r 2 )x(l) 



^(rOxCT) ^ /) (r 2 )x(t) 
VA(n)x(l) ^(r 2 )x(-L) 



Here <Pa/b( v ) describes an outgoing wave which has in 
principle the asymptotic form 



TA/B 



e ikr ^Y w (r A/B ); r AB = \r-R A/1 



with Yio(r A /g) denoting the spherical harmonic for I = 
1, m — 0, and k is given by Ti 2 k 2 /2 too = — 1 Ryd + huj. 
We have assumed linearly polarized Laser light with the 
quantization axis of Y\q(ja/b) coinciding with the axis of 
polarization. Moreover we have disregarded the residual 
charge left with each atom as part of the electronic charge 
is emitted. 

Although only the illuminated volume of atom A can 
contribute to the transition matrix element 

Mfi = 

/ / *}(ri,r 2 ) V perturb (r u r 2 ) %(ri,r 2 ) d 3 r 1 d 3 r 2 

J atom A J 



the final state ^/(ri,^) yields a current density 
j(r) = / [** f (r,r 2 )V*/(r,r 2 ) -c.c.]d 3 r 2 

h 
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J [*)(ri,r)W/(ri,r)-c.c] d 3 n = 
ja{v ~ Ra) + 3B{r - R B ) 



where Ja/b is associated with <Pa} b , and hence j(r) con- 
tains also a photo emission current coming from the non- 
illuminated atom B at a distance of 10 cm away from A. 
Similar considerations apply if one excites the molecule 
to a bound state which would spontaneously decay back 
then to the ground-state by emitting fluorescent light. 
If one repeats the excitation sufficiently often one would 
obtain as many fluorescence photons coming from the il- 
luminated atom as from the non-illuminated one. There 
is no experimental evidence that anything like that could 
ever happen. 

We are hence led to conclude that the concept of entan- 
glement (i. e. the Pauli exclusion principle when dealing 
with fermions) does not apply anymore if the atoms are 
separated by a macroscopic distance. The reason may 
be tracked down to the definition (J4j) of the probability 
density p(r) as the relative residence time that a particle 
spends in an elementary volume A 3 r around r, provided 
it is bound in a potential and thus occurs repeatedly in 
that volume. On pulling the two atoms of a H 2 -molecule 
gradually apart one arrives at a situation where one of 
the two electrons remains captured near the nucleus of 
atom A for a while, and accordingly the second electron 
stays captured near the nucleus of atom B for the same 
time. The Coulomb repulsion between the two electrons 
effects a correlated separation of the two electrons into 
the two regions^ If the inter-nuclear distance becomes 
large compared to the linear dimensions of the atoms, the 
time spans for tunneling of the "A-electron" (marked by 
the index "1") into the B-region and vice versa become 
enormously long compared to the time required to tra- 
verse the associated atom. The time T for the photo- 
excitation process will therefore be many orders of mag- 
nitude shorter than the tunneling time. Given this situ- 
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ation, the definition ([4]) yields 



1-2) 



PA/B{ r i/2) f° r r i/2 around nucleus A/B 
else 



where - different from Eg. (|133p - the densities Pa/b( t i/2) 
now integrate to unity. The two electrons do not appear 
entangled any more, and only the A-atom will now emit 
an electron under the exposure of light. 



XXII. DECOMPOSING AN EXPERIMENTAL 
SETUP INTO THE QUANTUM SYSTEM UNDER 
STUDY AND A REMAINDER. SCHRODINGER'S 

CAT 

One of the puzzling credos of the Copenhagen interpre- 
tation of quantum mechanics consists in the conviction 
that an experimental setup for performing measurements 
on microscopic particles, has to be subdivided "some- 
how" into the particles under study and a remainder 
that functions as a classical system. This decomposi- 
tion is known under the name "Heisenberg-cut" . Yet 
from an unbiased point of view it appears to be self- 
evident that an experimental setup as a whole represents 
a many-particle system each part of which is subjected 
to the same laws of quantum mechanics as the particu- 
lar portion that constitutes the object under study, an 
electron in a diffraction chamber, for example. We shall 
use this example to demonstrate the consistency of this 
standpoint, but we limit ourselves to considering a sys- 
tem that merely consists of just one specific apparatus 
plus a particle undergoing diffraction in it. The gen- 
eralization to the inclusion of the entire environment is 
obvious from the ensuing considerations. 
We assume that the system is made up of N particles, 
a subset consisting of atomic nuclei which we number 
by a label a, and N e electrons, one of which represent- 
ing the single particle of interest, the "test particle". To 
keep the notation simple, we limit ourselves to consid- 
ering only electrostatic particle interactions of the kind 
described by the many-body potential (|119p . If the test 
particle has left the cathode of the setup it is kept by elec- 
trodes, diaphragms and lenses at a macroscopic distance 
away from all kinds of surfaces it might strike and where 
it might get captured. Thus, the associated one-particle 
wave function tp e (f, t) which describes the electron on its 
way through the apparatus to the screen or detector, has 
de facto zero overlap with the wave function of the N — 1 
remaining particles of the apparatus. Still, in standard 
setups it is intended that the particle hits a secluded por- 
tion of material on its way to the monitoring device, a 
diffracting single crystalline foil of metal, for example. 
But in the majority of cases the contact time is so short 
compared to the electronic excitation times of the mate- 
rial that the test electron cannot mingle with the other 
electrons. Further below we shall briefly discuss promi- 
nent exceptions. 



Similar to the case of the H2-molecule with macroscop- 
ically distant nuclei, one is justified then in assuming a 
factorization of the total wave function 

*jv(r,r 2 , ...r N ,t) = 

ip e (r,t)V N -i(r 2 ,r 3 ,...r N ,t) (134) 

where the spin coordinates have been suppressed for sim- 
plicity. If we insert this into the A^-particle Schrodinger 
equation (|12ip we obtain 



^N-liUlpe + Ipeih'fy N -i = 



(JV-1) 



*AT_1 + *Ar_i HI 1p e + Vtotlpe *N-1 (135) 



where Vtot denotes the total (Coulombic) interaction po- 
tential between all particles 



in 



N 



V a pparatus(r2 + 



e 2 Z* 



J=2 



47r e r,- 



(136) 



and from Eq. fUg) 

N 
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+ VextXrj) 



N-l 



2 toq 



+ V ext .(r) 



=Hf, 



In Eq. (|136p \Zj\ stands for the number of elementary 
charges, i. e. 



Z< 



— Z a if j runs over the a*' 1 nucleus 
1 if j refers to an electron . 



If one multiplies Eq. (|135p by ^* N -i ano - performs an in- 
tegration with respect to r 2 , . . . one obtains 



V, 



apparatus 



•* 9 T 

ih— $ N -i- 
at 



d 3 r 2 . . . d 3 r N 



where V e (r, t) represents a one-electron potential defined 

as 



N 



V e (r,t) 
e 2 Z 4 



47r e |r 7 - 



^N-id 3 r 2 ...d 3 r N . (137) 



Since the bracketed expression under the integral in the 
above equation vanishes, we arrive at 

ih J^e(r) = [H e a +V e (r,t)]Mr) ■ 
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Thus, the wave function of the electron under study 
obeys, in fact, a one-particle Schrodinger equation. 
There are certain cases in which the contact time of 
the test particle is not short enough, and hence there 
is a non-vanishing probability that the particle mingles 
with those of the target. To get a rough picture of this 
situation, we describe the wave function instead of Q134p 
by 

%N{r,r 2 , ...r N ,t) = c (t)i/j e (r,t) fjv-iO^, r 3 , . ..r N ,t) 



d(t) ^ pt (r,r 2 ,...r N ,t) 



(138) 



where co(i) and c\(t) are real- valued functions with the 
property |cq(£)| 2 + | ci (*) | 2 = 1, in particular 



co(t) 



and hence 



|co(t)| 2 = 



IciWI 2 = 1 



(139) 



Here r refers to a characteristic interaction time with 
the target, and | Ci (*) | 2 is the probability with which the 
test electron is captured by the target. Thereby it loses 
its identity as the "test electron". The latter effect is 
expressed by the property of ^> c ^ pt being antisymmetric 
with respect to interchange of any two particles out of 
the set of N electrons. 

Inserting from Eq. (|138p into the Schrodinger equa- 
tion (|121[) we obtain 



c (t) *at_i ih tb e + il) e i%-$> N _ x 

ihto{t) ipe^N-i + c 1 (t)ihi>™ pt 

co(t) f 



+ ihc 1 {t) ^ pt 



c 1 (t)H^>'j apt 



■ N 



Vtotlpe &N-1 

(140) 



The functions ^n-i and ^ N satisfy the associated 
time-dependent Schrodinger equations 



H, 



(JV-l) 







V, 



and 



apparatus 



Vn-i (141) 



If we insert this into Eq. (|140p . multiply the result 
in front by '^* N _ 1 and perform an integration over 
r 2 , t*3, . . . rjv, we obtain 



co(t) 



H e - V e '(r,i)j Mr,t)+ 

ih r (/) / >l' v I'''-'- ■■■fx.!} 

c (r,r 2 ... r N ,t) d 3 r 2 d 3 r 3 . . . d 3 r N = . (142) 
Here we have used ih co(t) — — i co(t) and set 



x^ pt < 



m c oW = 

V e '(r,i) = V e (r,t) +iV e where 



Ve 

n 



i 

2t 



The imaginary part of V e '(r,i) is commonly referred to 
as "optical potential" . 

According to our classification of the electron 
under study as either "distinguishable" or "non- 
distinguishable" the associated total probability density 
p(r,t) splits (almost quantitatively) into the "either- 
and or-probability" 

p(r,t) = \c (t)\ 2 \Mr,t)\ 2 + 
v * ' 

=Po{r,t) 

\ Cl {t)\ 2 J \*% pt (r,r 2 ,...r N ,t)\ 2 d 3 r 2 ,d 3 r 3 ...d 3 r N 



=Pl(r,t) 



which means 



S e (r,t) = J *^_ 1 (r a ,r 3 ,...rj V ) 
x*™ pt (r,r 2 ,r 3 , . . . r N ) d 3 r 2 , d 3 r 3 . ..d 3 r N w 0Vr,t. 



It follows then from Eq. (1142[) that ip e (r,t) solves the 
modified Schrodinger equation 



(143) 



which describes situations one encounters, for example, 
in experiments on low energy electron diffraction (LEED) 
at surfaces of solids. 

The mod squared of the actually not completely vanish- 
ing overlap 

S(t) = f S e (r,t)re(r,t)d 3 r 



determines the transition probability 1/r. Eq. (|138|) and 
the resulting Eqs. (ll41| ) and (I143P are pivotal in describ- 
ing generic quantum mechanical processes, a subset of 
which plays the role of "measurements"—. For exam- 
ple, when one is dealing with a setup where an electron 
traverses the legendary double slit diaphragm, defined 
by the potential (|137[) , the function ^ c ^ pt (r, r 2 , . . . rw,t) 
describes the situation when the electron has been cap- 
tured by the detector which is a part of the "apparatus" . 
In spirit this in keeping with a statement by Hartlc and 
Gell-Mann^: "In a theory of the whole thing there can be 
no fundamental division into observer and observed" Our 
approach reflects even more directly the standpoint taken 
by v. Kampeni£: "The measuring act is fully described 
by the Schroedinger equation for object and apparatus to- 
gether... " 

It is the archetypal combination of a particular setup 
and "pointer readings" of a detector that enables the 
experimentalist to determine certain properties of the 
one-particle quantum system by extracting the sought- 
for information from the solution to the correspond- 
ing Schrodinger (or Pauli) equation that yields j(r) = 
p(r,t) v(r 7 t) at r detector- Eigenvalues of hermitian oper- 
ators can only be obtained via this detour, and for funda- 
mental, mostly experimental, reasons, only with limited 
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accuracy. 



and hence 



The paradoxical situation which one runs into if one en- 
dows the "observer" (or "measurer" ) with an unrealistic 
meaning, is illustrated by Schrodinger's cat example 47 : 
An alpha-particle emitted from some radioactive material 
triggers a device that kills a cat in a closed box by releas- 
ing a poisonous gas. Of course, the moment of radioactive 
decay does in no way depend on the particular setup. Our 
description of this process would be based on Eq. (|138p 
where ^Ar(r, ra, . . . r^, t) on the left-hand side now rep- 
resents the wave function ^ ga s+cat(t) of the system cat 
plus gas, ip e (r,t) has to be replaced by an iV-particle 
wave function ip po ison referring to the only weakly "cat- 
overlapping" molecules of the poisonous gas set free by 
the device, and 5 , Ar-i(r 2 , r 3 , . . . rj\r, t) has to be identified 
with the many-particle wave function ^ ca t of the live cat. 
After the elapse of a time ~ r the system's wave function 
^ gas+cat{t) has attained the form ^ C apture(t) where the 
poisonous molecules are now a part of the cat. It solves 
the time-dependent Schrodinger equation of the united 
system. The time-evolution of ^ C a P ture{t) describes all 
the atomic (chemical) processes that eventually lead to 
the cat's death. It is this time-dependent process that 
is familiar from ab initio calculations on chemical reac- 
tions. The latter are completely "self-controlled". There 
is definitely no "observer-induced" influence. From this 
point of view it appears to be rather absurd that ortho- 
dox quantum mechanics interprets Eq. (|138p with the ex- 
plained new meaning of the wave functions as a superpo- 
sition of a "live" and a "dead" -state of the cat, and only 
on opening the lid of the box by an observer, ^ gas +cat 
collapses onto the wave function of a live or dead cat. 



XXIII. THE ORIGIN OF PARTICLE SPIN 



We ignore here and in the following the minute departure 
of g from 2 due to quantum electrodynamical corrections. 
The radius of the rotating electron sphere was identified 
with the classical electron radius 2.8 • 10~ 13 cm. Lorentz 
immediately demonstrated to Uhlenbeck and Goudsmit 
that the electron mass would actually be larger than that 
of a proton if the magnetic moment of a Bohr magne- 
ton would be confined to that sphere. Moreover, the 
speed at the equator of the rotating sphere would by far 
exceed the velocity of light. Although these objections 
definitely disqualified the rotating sphere as a model of 
electron spin, it is still used, tacitly implied or appears 
concealed as "intrinsic property" in the analysis of most 
of the present-day experiments involving spin-orientation 
or spin flips. 

If an "eigen-rotation" cannot explain the occurrence of 
a mechanical spin moment associated with a gyratory 
electronic motion, what else can be responsible for it? 
The following considerations are based on the idea that 
"spin" is not a property of the particle but is rather 
a property of its quantum mechanical state. 
Our description of particle motion as modified by 
stochastic vacuum forces makes it particularly suggestive 
to correlate - similar to the explanation of zero-point mo- 
tion of oscillators - particle spin with the quivering mo- 
tion that results from those forces and vanishes as h tends 
to zero. (This applies, of course, to all the other quantum 
mechanical ground-state properties as well.) To illustrate 
this we consider the simplest case of a hydrogen electron 
exposed to a magnetic field B — Bz e z in its ground 
state ifjo( r )- In Fig [3] we show a schematic distribution 



In 1925 Uhlenbeck and Goudsmit^ suggested in a 
widely recognized paper that Pauli's idea^ of a fourth 
quantum number in the description of electronic states 
of atoms might be associated with the rotation of an 
electron about its own axis thus giving rise to an extra 
angular moment. From the analysis of atomic spectra it 
was clear that the magnetic moment generated by such 
a rotation of the electron as a charged sphere had to ex- 
actly equal the Bohr magneton 



■S 



fJ-B 



eh 
2mn 



There was also experimental evidence that the associated 
mechanical spin moment S - different from the atomic 
orbital momentum - would not obey the classical law of 
magneto-mechanical parallelism according to which \ib 
should differ from S by a factor 5^-. In actual fact this 



factor had been found to be 



\S\ 



instead so that 



h 
2 



FIG. 3: Spin effective components of the quivering motion 

of positions that the electron has successively taken at 
times ti and equal time intervals ti+i — U = At where 
At is very small compared to T. This time span has 
been introduced in Section IIIII in connection with defin- 
ing the probability density p{r). The z— axis is thought 
to run through the atomic center perpendicular to the 
plotting plane. At each of the points the electron pos- 
sesses a velocity which we decompose into a radial and 
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a z— component, and in a component perpendicular to 
the z— axis. Only the latter are indicated by arrows. For 
symmetry reasons there will be as many positive as nega- 
tive radial and z— components in the elementary volume 
around each point. They average out. We subdivide the 
set of arrows into two subsets associated with left-hand 
and right-hand circular motion, respectively. One might 
surmise that in the presence of the magnetic field one of 
the sets becomes empty and the other set now gives rise to 
a net circular current so as to minimize the total energy in 
that electronic state. As a defining property, the energy 
gain is proportional to B z and vanishes as h — > 0. As we 
know from Section IXIVI the magnetic field causes in the 
state ipo(r) a current density j(r) = — \ipo(r)\ 2 A(r), 
but the energy gain from this goes as B 2 because of 
AE = J j ■ A d 3 r and B = V x A. The omission of 
the empty subset of arrows does not at all change the 
distribution of points defining ipo(r). The energy gain is 
provided by the vacuum fluctuations in complete analogy 
to the zero-point energy of oscillators. 
We thus arrive at the conclusion that the circular cur- 
rent which occurs on allowing the quivering motion of 
the electron to become asymmetric does not change the 
probability density which is characteristic of real-valued 
solutions to the Schrodinger equation. But it definitely 
yields a physical effect that has so far been outside our 
formal framework. In discussing certain properties of so- 
lutions to the Dirac equation Schrodinger— was led to 
a similar interpretation of particle spin and coined the 
irregular particle motion causing it "Zitterbewegung" . 
The additional spin-dependent interaction with a mag- 
netic field occurs also in complex-valued states ip(r) — 
\ip(r)\e lip ( r \ In that case there is an additional set of ar- 
rows superimposed on those shown in Fig.l. That set 
consists of arrows depicting v(r) — ^ Vip(r) at the 
various points distributed according to |t/>(t")| 2 . Clearly, 
a linear superposition of those arrows is only possible 
as long as the velocities are within the non-relativistic 
regime. Otherwise the superposition is affected by spin- 
orbit coupling as a result of which tp^ (r) and (fi (r) now 
become different. This point will be taken up again in 
Section KxrXl 



XXIV. GENERALIZING ONE-PARTICLE 
QUANTUM MECHANICS BY INCLUDING 
PARTICLE SPIN 



Several suggestions have already been made to incor- 
porate particle spin into a theory that is akin to the ideas 
of the present article (s. e. g. Dankel^, Dohrn et al^, 
Nelson—). We believe, however, that our approach of- 
fers - in the spirit of a statement by v. Weizsacke r 52 i 82 - 
"something immediately comprehensible" . 

The points associated with the two subset of arrows in 
Fig.l define probability densities p-[(r) and Pi(r) with 
| and J, referring to the respective direction of the spin 



moment. Both densities integrate to unity 

p m) (r)d 3 r=l. (145) 



To keep the formalism flexible at the outset we consider 
a situation where the total probability density is not yet 
a pure "up" or "down" density 



p(r) = \a\ 2 p^r) + \b\ 2 p l (r) 



(146) 



with a and b denoting coefficients whose modulus squares 
sum up to unity 



\a\ 2 + \b\ 2 = l 



(147) 



It is obviously not possible to partition the wave function 
analogously: ip(r) = aip^(r)+btpi(r) because ip*(r)ij}(r) 
would contain cross-terms. However, if one introduces a 
two-component spinor of the form 



ip(r) 



oV'tW (o)+^lW (JW) 



and its adjoint ^ (r) = (a*ijj*(r), b*ip*(r)j where 

\^ m (r)\ 2 d 3 r = l, (149) 

one obtains as intended 

f (r)^(r) = M 2 \^(r)\ 2 + \b\ 2 |^(r)| 2 = p{r) 
and 

^(r)jP(r)d 3 r = 1. (150) 

We consider the Bohr magneton as known from the ex- 
periment. Thus, the energy densities of the interaction 
with the magnetic field for "up"- and "down" -spin may 
be cast as 

-p B B z \a\ 2 ^(r)^(r) and +pbB z \b\ 2 ip*(r) ipi(r) . 
from which the total interaction density results as 

U,nagn. (r) = (r) p B R ^_{r) (151) 

where we have introduced a matrix 
B z 



B 



-B, 



(152) 



Likewise, we may cast the non spin-dependent energy 
density of the electron as 



where 



H = H Q + V(r) and H Q = ^ ^ A (r)f ^ 

2 rrio 

and with V(r) denoting some potential in which the elec- 
tron moves. 
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XXV. THE TIME-DEPENDENT 
NON-RELATIVISTIC PAULI EQUATION 

The basic two constituents of our approach, 
viz. \ip(r,t)\ 2 and Vy(r, t) remain unaffected by 
our incorporation of spin. Hence, it is completely in line 
with the conceptual idea of our approach to assume that 
the two theorems of Ehrenfest stay unaffected as well. 
That means, according to Ehrenfcst's Second Theorem 

(■w) = ^- ) = y il{r,t)riAr,t) 

+^(r,t)r<if)(r,t) d 3 r, (154) 
and we have alternatively from Eq. (|94|) 

(v) = — f ^(r,t)PiP(r,t), (155) 
mo J - - 

which holds without modification also for the spinors we 
have introduced. Exploiting the relation 



[H r-rH ]^(r,t) 



-i — Pip(r,t), 
m — 



which follows from simply applying the chain rule, we 
may combine Eqs. (|154p and (|155[) to obtain 



/ 



Ho + ih 



ij}'(r, t) r 



d_ 

dt 



ip\r, t)^j ri/j(r, t) d 3 r - 
ip{r,t) d 3 r = Q. 



Hn — ih 



dt 



(156) 



This equation holds for any t if tp(r, t) satisfies 



4(r,t) = -D4>(r)-F(r,t)il>(r,t), 



(157) 



and correspondingly 
d 



H Q + th- 
at 



tp^(r,t) = -ift (r)D-ip*(r)F(r,t) 



where F(r,t) is some integrable real-valued function 
and D_ denotes some unitary 2x2-matrix that will be 
specified later to meet requirements of Ehrenfest 's first 
theorem. 

The expectation value of the force exercised on an 
electron which moves in a potential V(r) and simulta- 
neously - through its magnetic moment - feels a force in 
a spatially varying magnetic field B(z,t) — B z (z,t) e z 
may be cast as 

(F) = -J^(r,t) [V{V(r) + p B R}] i L (r 7 t)d 3 r 



(eA(r,t)) 



(158) 



The appearance of the induction-derived force 
— {eA(r,t)) is a consequence of Eq. (|9"T]) . 



We perform an integration by parts on the first integral 
and obtain 

(F) = J [Vj&^r.t)] {V(r)+p B R}i L {r,t)d 3 r+ 
J ^{r,t){V{r) + fi B B} V^(r,i) d 3 r 

-(eA(r,t)). (159) 

From Eq. ([94| we have 



\p) = — / ^ f (r, t) [-ihV - e A(r, t)] ip(r, t) d 3 r 



which we rewrite 
(P> = 



Vif)d 6 r+ I 7/i f V 



-ih^- ib ) d 3 r 
dt- 1 



On forming (F) 
we obtain 



=-/VV.t(-iK^) d s r 

ip^ if; eAd 3 r . 
(Ehrenfest 's First Theorem) 



/ (V^ f ) {V + p B R- th§- t } ^d 3 r+ 

J (+ih^^V^ + £{V + fi B B}V-iP d 3 r = 0. 

(160) 

If we here eliminate the time-derivatives using Eqs. (|157|) 
and identify D with p B R this equation takes the form: 

- j\iV^)Ho±+{H Q £)V^\ d 3 r+ 

J (V^ f (r,t))?A(r,i) +^ f (r,i) V±{r,t) 

x{V(r)- F(r,t)}d 3 r = 0. 
Since the first integral vanishes we arrive at 

Jvp(r,t){V(r)-F(r,t)}d 3 r = V*. (161) 

We first consider the possibility that the expression in 
curly brackets does not vanish, but the integral does. 
As we have emphasized in defining probability densities 
p(r,t) and average velocities v(r,t) through Eqs.Q and 
Q , non-stationary states require a certain sample-time T 
of the particle under study to allow its time-derived prob- 
ability density p(r, t) to become quasi-stationary. Hence, 
if we introduce at t — to a small perturbational potential 
V(r) -> V(r) + Sv(r,t) where t < t < T, the proba- 
bility density p(r,t), and thus its gradient remain prac- 
tically unaffected, but the bracketed expression is now 
definitely different. We are hence led to conclude that 
Eq. (|16ip can only be satisfied if F(r, t) = V(r) holds for 
any time. That means - because of Eq. (|157p - that the 
spinor function ip(r,t) solves 

H + V(r) + p B Bj V>(r, t) = * ^ ±(r, t) (162) 
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with Hq as defined in Eq. (|153[) . This constitutes the 
time-dependent non-relativistic Pauli equation. 



XXVI. THE CAYLEY-KLEIN PARAMETERS 
AND PAULI SPIN MATRICES 



We want to adapt Eq. (|162|) to a situation where the di- 
rection of the magnetic field no longer coincides with the 
z-axis of the coordinate system. This can be achieved by 
exploiting a surprising alternative to the standard form 
of rotating the coordinate system by applying orthog- 
onal 3x3 matrices. The idea goes back to Felix Klein 
(S. Goldstein^) and is related to earlier work of Cay- 
ley. He considers the rotation of the coordinate system 
(x,y,z — > x',y',z') to be performed in three steps de- 
scribed by the Euler angles <fi, 6 and ip shown in FigJU 
Instead of representing the position vector r by a col- 






FIG. 4: Euler angles 



umn matrix he uses a 2x2-matrix P(r) of the form 



z x — ly 
x + iy —z 



(163) 



In place of the standard 3x3-rotation matrix one now 
has a 2 x 2-unimodular matrix 



whose elements - the so-called Cayley-Klein parameters 
- arc connected to the Euler angles through 



= ie^ sin 



7 
8 



ie -W-*) sin . 



e -i(*+*) cos : 



(165) 



After the three steps of the rotation have been performed 
the original position vector r = (a;, y, z) is now associated 
with the new coordinates x' , y' , z' that may be obtained 
from the transform 



QEQ = P (x ,y ,z') = 



x — i y 
-z' 



(166) 



x + iy 

where Q + denotes the adjoint of Q, and we have 

£+£ = ££+ =1 (167) 
The matrix B_ had been defined in Eq. (|152|) as 
B z \ „ / 1 



= -\ -B z 



= B Z 



-1 



(168) 



In Klein's representation the point r = (0, 0, z) attains 
the analogous form 



P(r) = z 



1 

-1 



Hence, B_ has to be required to transform under coordi- 
nate rotation as P: 



R' = QBQ + 



(169) 



For a general orientation of the coordinate system with 
respect to the magnetic field B_ has the form analogous 
to P in Eq.pSHJl, viz. — 



B 



B z 



B x -iB, l 



B x + i By —B z 
This matrix can be decomposed 



B = B r a + B„ a 



B z <7 



(170) 



(171) 



where the three matrices on the right-hand side are just 
the Pauli spin matrices 



1 

1 



a = 
= y 



-i 

1 



1 

-1 



(172) 



They are commonly lumped together in the form of a 
vector 



a = a e T + a e v + a e. 



(173) 
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The matrix B_ in Eq. (] 1 Tl [) may therefore be cast as 

11=3 -B. (174) 
The Pauli equation (|162[) then attains the familiar form 

Ho + V(r) + ng a ■ b] ±{r, t) = ih^ ±{r, t) . (175) 

Actually, the spinor in this equation should be marked by 
a prime because it has changed under the transform as 
well. We have dropped the prime for simplicity. Since the 
density of the magnetic interaction energy is, of course, 
invariant under rotation of the coordinate system 

u magn.{ r ) = u magn.( r ) > 

it can be shown then that the new ip' is connected to the 
original ip through 

= Q± (176) 

and correspondingly 

This becomes obvious from forming 

figf (= : u' mag (r',t)) = ^q+gqi. 

If we insert Eq. (|169p on the right-hand side we obtain 
^Q + gQiP = ^ Q + Q BQ + Q ip = u mag (r, t) . 



Because of 



the probability density is also invariant under rotation of 
the coordinate system which is consistent with our idea 
of a spin-defining motional decomposition at the begin- 
ning of our considerations. Moreover, if the state of the 
particle in the original coordinate system has the form 



or ipi(r) 



it becomes after coordinate rotation 
^ T (r)=^ t (r)£(r) (l 



(177) 



a(r) 



/3* (r) 



f x (r) = ^(r)£(r) 




1 



p{r) 



+ a*(r) 



where a(r) and (3(r) are the Cayley-Klein parameters 
describing the rotation which we have allowed here to be 
different at different positions r. 

The spin orientation with respect to the direction of a 
magnetic field is already uniquely defined by the two an- 
gles 9 and 4>. Hence one is at liberty to choose -0 at 
will without loss of generality It is convenient to set 
ip = — 7r/2. We consider the projection of the unit vector 
e' z onto the original x/y-plane where it makes an angle 
tp with the x-axis. This angle and the Euler-angle <f> are 
interrelated 



^ + 2 



If one inserts this relation into Eqs. (ll64p and (jl65[) . Q 
takes the familiar form 



_ / cxp[§^(r)] cos^ exp[-|ip(r)] sin^ \ 

— ^ \ N l M • e ( r ) r i I M SM I ' 

— \ — exp[^ip{r)\ sin-j 1 exp[—^ip(r)\ cos / 

(178) 

The above considerations on the magnetic interaction en- 
ergy starting with the expression (|151|> carry over to the 
spin momentum 



6(r) 



Da| 2 |^(r)| 2 -|6| 2 |^(r)| 2 ] d 3 r. (179) 



The symbol S z refers to the effective spin moment in the 
z-direction with respect to which the functions ipm)( r ) 
have been defined. In complete analogy to (|168p this 
expression can be compactified by introducing 



so that 



= 2 { -1 



(S z ) = / ^ + (r)S *P(r)d 3 r. 



(180) 



(181) 



In case that the functions tp^n-\(r) refer to a z'-direction 
that belongs to a rotated coordinate sytem x',y',z', we 
have in analogy to Eq. (|169p 

g z , = £^£ + - 

If we use the analogous relations pertaining to Eqs. (|170p 
up to (|174p we may cast S , as 



S 



- [a x a + a y a + a z a 
2 — x — y — ; 



(182) 



with a x , a y , a z denoting the component of the unit vector 
e z r in the z'-direction 

a x = cos ip sin 9 
a v = sin if sin 9 
a z = cos 9 . 
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It is convenient to introduce a vector S (commonly re- 
ferred to as "spin operator"), which is analogous to a, by 
setting 



S=-a. 



(183) 



Eq. (|182[) may then be cast 

= z ' = ez ' ' ' 

and hence we have 

(5,/) = 1 ^ + (r) £ z , V;(r) d 3 r = e z , ■ (S) 

where 

(S) = [ i/j + {r)STp(r)d 3 r. (184) 



If V>(r,t) has the form ([177]) . Eq.(TH3]) yields (5) = 
±| e z . On the other hand, if ip(r,t) possesses two non- 
vanishing components, there will always be a coordinate 
system that is rotated with respect to the present one, 
in which (S) becomes ±| e z >. One only has to turn the 
pertinent z'-axis in the plane spanned by the original di- 
rection of (S) and the original z-axis until e z i is parallel 
or anti-parallel to (S). 



XXVII. SPIN PRECESSION IN A MAGNETIC 
FIELD 

So far we have assumed the magnetic field and the 
spin direction to be collinear. As an example for a non- 
collinear situation we consider an electron that is bound 
within an atom where it is initially exposed to a mag- 
netic field along some direction. We omit here discussing 
the details of its spin alignment due to some minute time 
dependent perturbations and simply assume that it has 
eventually attained a stationary spinor state in which 
its spin momentum points parallel or anti-parallel to the 
direction of the magnetic field. If one now changes non- 
adiabatically the direction (and in general inevitably also 
the magnitude) of the magnetic field, the spin momentum 
can - without an appropriate external torque - not adjust 
to the new field direction, and hence the previously ex- 
isting collinearity no longer obtains. As we shall show by 
discussing the pertinent solution to the time-dependent 
Pauli equation (|175p . the spin momentum now precesses 
about the new direction of the magnetic field in a com- 
pletely classical way. 

We identify the initial direction of the magnetic field with 
the z'-axis of a "primed" coordinate system in which the 
spin-aligned state of the electron has the form 



^'(r')=^(0 



(185) 



where ip' (r') is the energetically lowest lying solution 
to the Schrodinger equation of the one-particle system 
under study. We denote this solution by ipo(r) in the 
unprimed coordinate system in which the new magnetic 
field lies along the z-direction and in which the spinor 
(|185|) can be cast as 



' <ip'(r') = ijj (r) e l ^ cos- 



—ib (r) 



(186) 



—lb (r) 



where 9, <p, ip(= 0) are the Euler angles that refer to the 
interrelation (#', y' , z') — » (x, y, z)). Hence we have 



±(r) =^(r)+'tp JOl (r). 



(187) 



Note that the unit spinors in Eq. (|186p are now referenced 
to the new z-axis! 

We now consider the Pauli equation (|175[) for the time- 
independent case in the absence of a magnetic field in 
which case ^^(p) are independent degenerate solu- 
tions and ipo(r) satisfies the associated Schrodinger equa- 
tion 

H(r)ipa(r) = E ipo(r) . 

For B z 7^ the two spinors belong to different energies 
-E/otC!) -fti fJ-B B z and their sum does not satisfy the 
time-independent Pauli equation any more. However 

^(r,t) = V^ OT (r)e"^ £;o T t +^ oi (r)e-^ Eo i t (188) 



solves the time-dependent Pauli equation (|162p if we dis- 
regard effects of second and higher order in the magnetic 
field. We now insert the definitions of i/j^^^r) from 

above and obtain 
4>(r,t) = ifj Q (r 



COS - 
2 



dl . 6 fO 
sin - 
2 U 



e~i Eot . 



(189) 



Here we have made use of -Eotd) = ± Ms B z and 
introduced the frequency ujl which is defined through 



Sot - Eqi = 2 [ibB z = Tilol ■ 
In complete analogy to Eq. (|186p we form 

and calculate the expectation value of S 
(S) = J iP + (r,t)Si>(r,t)d 3 r = 



(190) 
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= jy + (r',t)gSg + i/{r',t)d 3 r' 

= | [C0s((y3 — UJLt) Sin0ea; + 

sin(y) — UJLt) sin e y + cos 9 e z 



(191) 



Thus, the vector (S) of the spin momentum moves on a 
circular cone with an apex angle of 2 9 about the direction 
of the magnetic field and its projection onto the x/y- 
plane rotates at an angular frequency ujl, the "Larmor 
frequency", about the z-axis. According to Eq. (|190p this 
frequency is given by 



HbB z 

n/2 



(192) 



The spin precession is completely analogous to that of a 
classical spinning top which rotates about its symmetry 
axis at an angular frequency ui and is exposed to the 
gravitational field of the earth. The precession frequency 
lup is in this case given by 

Fr s 
UJP = — ' 

where L denotes the absolute value of the angular 
momentum, F is the absolute value of the gravitational 
force acting on the top's centroid, and r s is the distance 
of the centroid from the point of support. In case one 
has instead of a gravitational field a magnetic field and 
if the spinning top possesses a magnetic moment fip, 
one has Fr s = fisB z . Inserting this into the classical 
equation for u>p and setting L — h/2 one obtains exactly 
the expression (|192p for the Larmor frequency. If the 
magneto-mechanical parallelism would also hold for the 
spin momentum, that is if g in Eq. p44[ ) were equal 
to one, the magnetic moment would be ^/is, and the 
precession frequency would be smaller by a factor 2 in 
striking disagreement with the experiment. 
The completely classical behavior of a precessing spin 
moment in a magnetic field can also be made evident by 
the following consideration. 

Using Eqs. (|144p and (|19ip we may express the spin- 
derived magnetic moment Ms P i n as 

Mspin = Ms [cos(v? - UJLt) sin0e x + 

sin(<p — UJLt) sin 9 e y + cos 9 e z ] . 

The time derivative of Eq. (|191|) can be written 

j| < S >= [i B B z [sin(<ys - uj L t) sin6»e x ™ 

cos(ip — uiLt) sm6e y ] , 



where we have used ujl 
B = B z e z and 



(193) 

2fiB z /h. We observe that 



Hence, the right-hand side of Eq. (|193[) can be cast as 
/is B z [sin(<p — UJLt) sin^e^ — cos(<^ — uiLt) sm6 e y ] = 

M S pin x B . 

The result may be written 



dt 



< S >= Mspin x B . 



(194) 



This is identical with the classical equation of motion de- 
scribing the temporal behavior of a spinning top that is 
acted upon by a torque Mspin x B. It corresponds to 
Ehrcnfest's First Theorem, and it is this equation (|194[) 
which governs the phenomena encountered in electron 
and nuclear spin resonance. (S. e. g. Slichter— .) In ap- 
plying magnetic resonance techniques one has to sup- 
plement Eq. (|194[ ) by perturbational terms that cause a 
change of the precession cone. An equation of this kind 
was put forward by Bloch 58 in 1945. If the atom is not 
exposed to a time-dependent perturbation the spin keeps 
precessing on the cone without changing its apex angle 
even when the strength of the magnetic field adiabati- 
cally increases or decreases. A change of the absolute 
value of B only changes the Larmor frequency ujl- 

As opposed to the impression that is commonly invited 
by even the most recent literature, Eq. (|194p constitutes 
a purely quantum mechanical result and is in no ways 
"semi-classical" or "macroscopical" . The fact that from 
our derivation |(S%| may attain any value, seems to 
contradict the principle of "orientation quantization" ac- 
cording to which | (5% | may equal only integer multiples 
of h/2. Clearly, if (S) z is not parallel or anti-parallel to 
B but rather precesses about the direction of the latter, 
the electron emits magnetic dipole radiation until its spin 
is aligned. But this is a weak interaction, and therefore 
the state of non-alignment may well be regarded as meta- 
stable in certain experimental situations. 

Spin precession in a magnetic field exhibits a peculiar 
feature that relates to the occurrence of the argument 
^ in the exponential functions of Eq. (| 1 86[) . To see that 
we assume tp( r ) 1° represent a wavepacket of a free par- 
ticle that traverses a homogeneous magnetic field in an 
orthogonal direction. When the wavepacket enters the 
magnetic field the spin component perpendicular to the 
field may point in the x-direction which is also the direc- 
tion of flight. We then have 

ip = and hence e ±l 2 = 1 . 

During the flight 9 stays constant. When the wavepacket 
leaves the magnetic field after a full precession period we 
have 

(p = 2tt which means e ±l % = —1 . 

Hence 4'( r ) nas changed its sign, or one may just as well 
say, its phase has been shifted by tt. However, as can 
be seen from Eq. (|19ip . (S) points in the same direction 
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as at the beginning of the precession. This phase shift 
is well detectable in double-beam experiments with spin- 
polarized neutrons (s. e. g. Rauch^, Werner et al.— ). 



XXVIII. A THEORY OF THE 
STERN-GERLACH EXPERIMENT 

"...Phenomena of this kind made physicists despair of 
finding any consistent space-time picture of what goes on 
the atomic and subatomic scale... many came to hold not 
only that it is difficult to find a coherent picture but that 
it is wrong to look for one..." 

John Bell 61 

Quite a few attempts have already been made on a the- 
ory of the Stern-Gerlach (SG-) experiment^ 2 -. For a re- 
cent rather complete update of the pertinent literature 
see Home et al.— . But a coherent picture of the funda- 
mental mechanism is still missing. Most physicists seem 
to favor the idea that the electronic state of the atom on 
entering the magnet constitutes a linear combination of 
spin states "up" and "down", and the modulus square of 
the associated coefficients defines the probability of the 
atom for being either pulled up or down, that is par- 
allel or anti-parallel to the magnetic field gradient. On 
detection of the atom in the "up"- or "down" -beam the 
atomic wave function collapses onto the respective com- 
ponent of the linear combination. From our point of view 
this is unjustifiably associating the process of detection 
with some mystical influence of "observing", based on 
pure claim: the atomic beam would behave differently if 
it would not be detected. By contrast, we believe that 
the outcome of the experiment is completely determined 
by the time-dependent Pauli equation and is hence a re- 
sult of a "quantum mechanics without observer" . 
Our approach implies a linear combination of spin states 
as well, that is we describe the electronic ls-state of the 
atom that we shall consider below by 



a l 



where v denotes the velocity of the atom, and the coeffi- 
cients a-|- , a i have the property | a j | 2 + 1 ajj 2 = f. The unit 
spinors are referenced to the direction of the field gradient 
e z . Hence, the expectation value of the force acting 



3B_- 
dz 



on the atom in the SG-magnet is given by Eq. (|158p if we 
neglect the induction derived term and assume electro- 
static forces being absent 

(F atom ) = MB J tl to Jr, t) ^ ± atom {r, t) d\ e z . 

For simplicity we equate the field gradient to a constant 
so that (F a tom) reduces to 

QB 

(F atom ) = MB ^-i[|a T | 2 -H 2 ] e z . 



It can obviously attain any value between —(1b ^§r e z 
and +fiB e z depending on the value of the coeffi- 
cients when the atom enters the magnet. Therefore a 
splitting into two well separated beams cannot possibly 
occur as long as there is no particular mechanism which 
inhibits a random distribution. In the following we shall 
outline such a possible mechanism. 

We assume that the reader is sufficiently familiar with the 
essential features of the experimental setup. To simplify 
the line of argument we content ourselves with consider- 
ing the experiment by Wrede^ who used a primary beam 
of hydrogen atoms in a setup that was practically iden- 
tical with that of Stern and Gcrlach. Hydrogen offers 
the advantage of reducing the spin-orientation problem 
to that of a single electron. The standpoint we take here 
is akin to that of Mott and Massey-^ who remark: "From 
these arguments we must conclude that it is meaningless 
to assign to the free electron a magnetic moment. It is 
a property of the electron that when it is bound in an S 
state in an atom, the atom has a magnetic moment. "— 
The hydrogen atoms effuse from some source where they 
are (almost unavoidably) exposed to the terrestrial mag- 
netic field or at least to the weak fringe field of the SG- 
magnet. That field causes a weak Zeeman-splitting of 
the spin up and spin down level of the electronic ls-state. 
Because of the weakness of the splitting the two Zeeman- 
levels are at the temperature of the source equally occu- 
pied, that is, 50% of the effusing atoms have their elec- 
tronic spins oriented parallel to the weak external field, 
the spins of the remaining 50% atoms are anti-parallel. 
As the atoms approach the SG-magnet they feel in a co- 
moving coordinate system a magnetic field whose field 
strength increases continuously and will in general change 
its direction. We assume for simplicity that the spin ori- 
entation is transverse and that the atom moves along the 
x-axis of a laboratory-fixed coordinate system so that 
changes of the spin orientation will only take place in 
the y/z-plane parallel to the respective plane of the co- 
moving coordinate system. As soon as the field direc- 
tion in the co-moving coordinate system departs by a 
small angle 59 from the original direction of B = B z e z at 
the onset of the atom's trajectory, a small y-component 
By = B z sin 59 e y of the field appears as a consequence 
of which the magnetic moment of the atom experiences a 
torque —/is B z sm 56 e^,, where e v denotes the unit vec- 
tor in the direction of increasing azimuth angle ip in the 
x/y-plane. This torque causes a change L of the spin 
angular momentum 



- — sin( 
2 



'Wi e. 



where we have used 2[IbB z = Tiujl (Eq. (|190[> ). Hence, 
the spin momentum starts precessing about the new di- 
rection of the magnetic field. We ignore the slight tilt 
of the co- moving new cc/y-plane perpendicular the new 
field direction. 

We envisage a short time span for which we assume the 
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changes of 9 to be small so that 



sine 



S9 = 9t. 



(195) 



where t = coincides with the beginning of the rotation 
of the field. The following considerations exploit the typ- 
ical experimental condition that the precession frequency 
u>l is some orders of magnitude larger than the speed of 
the field rotation. (In the terrestrial magnetic field of 
magnitude ~ 5-10 T the precession frequency of the 
electronic spin is about 10 6 s -1 . At an atomic speed of 
10 cms , a distance of about 10 cm and a maximum ro- 
tation angle of tt/2 one has 9 «10 s .) As will become 
apparent from the following calculations we may limit 
ourselves to a short time span comprising only few pre- 
cession periods during which the magnetic field rotates 
only by a small angle (9 <§; 2tt) so that one is justified in 
assuming 9 to be constant: 



9 = const. 

The unit vector e v may be decomposed 
e v = —e x sin ip + e y cos ip . 



(196) 



At t = we have <p[t = 0) = — 5, that is e v — e x . Thus, 
it is advisable to replace ip with ip + -jj, but we omit 
denoting the new azimuth angle differently. Hence we 
have <p = for t = 0, and we obtain instead of Eq. (|196p 

e<p = e x cos ip + e y sin ip . 

The spin precession that now occurs is anti-clockwise 

ip = —lul that is ip — —lol t ■ 

Thus 

L = — 9 u>l \e x t cos u>L,t — e v t sin ul t] ■ 
2 

This results in a change of the angular momentum after 
one precession period T = 2tt/u>l 



Ujj B x Jq t COS OJl t dt— 



i ujl Jq t sitiljl t dt 



Hence, using 



27r p2ir 

£ sin£ d£ = — 2tt and / £ cos£ d£ = , 
ii Jo 



we may AL(T) cast as 

al(t) = \er^e v 

that is in the spirit of our approximation (|195[) 



AL(T) = - sin A9 e y . 



The y-component of the magnetic field which equaled 
zero at the beginning of the rotation is now given by 
By = B z sinA^ej,. That means: after one precession 
period T the magnetic field and the atomic spin angu- 
lar momentum have turned by the same angle A9. The 
spin orientation follows the magnetic field - within the 
present approximation - without slip, that is adiabati- 
cally. (This is similar to the physics of a spinning artillery 
shell whose spin axis follows the course of the shell's 
bending trajectory leaving only a small precession angle.) 
Thus, the atoms enter the SG-magnet (almost) fully ori- 
ented with respect to the SG-magnetic field. This applies 
to the atoms with anti-parallel spin orientation accord- 
ingly. Hence, the two beams leaving the SG-magnet re- 
flect merely the two kinds of atoms associated with the 
two Zeeman levels before they leave the reservoir. 
It is worth mentioning that Lev>^ carried out Stern- 
Gerlach-type experiments using beams of Na-, K-, Zn-, 
Cd- and Tl-atoms instead of Ag-atoms. The Zn- and 
Cd-atoms possess two s-valence electrons which results 
in a zero net spin momentum of the atoms and conse- 
quently one does not observe a beam splitting in the 
Stern-Gerlach magnet. On the other hand, Tl-atoms pos- 
sess a 6p-valence electron that is subjected to spin-orbit 
coupling. This gives rise to a Lande factor g — | as 
a result of which the effective magnetic moment is for 
Mj = \ given by 

Mf = VBg Mj = ^ hb ■ 

This is, in fact confirmed by the experiments. 
If one were dealing with atoms that possess a total angu- 
lar momentum J = (I ± |) % associated with 21 + 2 dif- 
ferent magnetic quantum numbers Mj, one would have 
21 + 2 different states in the initial weak field and there- 
fore as many different sorts of atoms entering the Stern- 
Gerlach magnet where they are deflected according to 
their magmetic moment. That means one would have 
21 + 2 different beams instead of 2. 

Our explanation of the SG-experiment is much in the 
spirit of Stern's conjecture that the spin of an atom re- 
sponds adiabatically to the directional change of the mag- 
netic field in which it has originally been aligned. In co- 
operation with Phipps^i he devised an experiment where 
one of the beams at the exit of a first SG-magnet was fo- 
cused into a linear set of three successive magnets whose 
weaker, essentially homogeneous fields pointed in three 
different directions perpendicular to the atomic trajec- 
tory. The difference between these directions was 120°. 
If the spin of the selected beam was pointing up after 
leaving the first SG-magnet and assuming that the spin 
would adiabatically adjust to the local magnetic field on 
its passage through the three magnets, it was thus to 
be expected that it would be finally back to its previous 
"up"-orientation. To test this the beam was sent into a 
second SG-magnet identically oriented as the first. There 
was only one beam coming out of this magnet indicating 
that the spin was pointing again in the same direction as 
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on entering the three "turn magnets". In other words: 
even after a turn of 360° no slip between spin orienta- 
tion and the direction of the magnetic field had occurred. 
We mention here only in passing that our result on the 
Phipps-Stern experiment agrees with that of Rosen and 
Zener— published already in 1932. Different from our 
more summary analysis these authors attempt to stay 
close to explicitly solving the time-dependent Pauli equa- 
tion. 

Surprisingly, the interpretation of the SG-experiment as 
demonstrating a coherent splitting of the de Broglie-wave 
of the incoming atom into two beams has become the 
most popular view on which a host of considerations 
on "measurement" is based. Papers on the so-called 
"Humpty-Dumpty-problem" (s. e. g. Englert et al. 69 ) 
deal explicitly with a possible reconstruction of the orig- 
inal single wave by appropriately merging the two coher- 
ent beams at a spot reached later. We believe that such 
thought experiments are without substance. As we have 
clearly demonstrated, the SG-magnet does not cause a 
splitting of the incoming matter wave. The SG-situation 
is distinctly different from that in neutron spin-flip ex- 
periments by Rauch and coworkers^ where a transverse 
spin polarized beam of neutrons hits a plate of a Si single 
crystal such that each matter wave packet splits up into 
two widely separated beams of packets due to dynamical 
diffraction within the crystal. This diffraction process is 
spin-independent. The two beams are coherently merged 
then by dynamical diffraction at a second Si-plate. 

Many authors give the impression as if there were not a 
shadow of doubt that Stern-Gerlach experiments with 
charged free particles (like electrons) are just as fea- 
sible as with spin-carrying neutral atoms. Bohr had 
very early pointed out (s. Wheeler and Zurek^i) that 
such experiments could not possibly succeed because 
"the Lorentz force would inevitably blur any Stern- 
Gerlach pattern" . Nevertheless, the literature on EPRB- 
(Einstein-Podolsky-Rosen-Bohm) correlation with pairs 
of fermions in a singlet state (s. e. g. Einstein et al4£, 
Bohm 44 ) abounds with allusions to "measuring sepa- 
rately the x/y/z-spm components" of the particles by 
means of Stern-Gerlach magnets. (S. e. g. Wigner— .) 
Even when one were dealing with neutral fermions what 
kind of mechanism should yield such information on those 
spin components? How would the time evolution of the 
respective solution to the time-dependent Pauli equation 
look like in this case? 



XXIX. THE TIME-DEPENDENT DIRAC 
EQUATION 

In trying to extend the theory to relativistic systems we 
retain the following two fundamental assumptions that 
characterize the non-relativistic quantum mechanics we 
have been dealing with so far: 

1. The universal existence of stochastic forces that ne- 



cessitate an ensemble description of the one-particle 
system under study. The fundamental constituents 
of this approach are: p(r,t) for the occurrence of 
the particle at r and time t and p(r, t) for the asso- 
ciated ensemble average of the particle momentum 

2. Lumping together the two real-valued functions 
p(r,t) und p(r,t) in the form of a complex-valued 
function ip(r,t) 

V»(r, t) = Vp(tv0 e l * {r ' t] (197) 

where 

p(r,t) =hVtp(r,t) . (198) 

From ip(r,t) = \ip(r, t)\ e 1 ^*"'*) one then obtains 
the momentum current density 

j p {r,t) = p(r,t)p(r,t) = 

±[iP*(r,t)prp(r,t)-iP(r,t)p ^ (r , *)] 

(199) 

where 

p = — i fi, V. 

Eq. (|198[) implies that p(r,t) is curl-free, that is, the 
stochastic forces do not cause friction. 
From Eq. (|199p follows for the expectation value of the 
particle momentum 

<p(t)>= J j P (r,t)d 3 r = J ^*(r,t)p4>(r,t)d 3 i200) 

If one replaces ip(r,t) with its Fourier integral 

^(r,*) = (27r)-f J C(k,t)e ik r d 3 k, 

one obtains on insertion in Eq. (|200p 
< p(t) >= J tJj* (r, t) p ip(r, t) d 3 r = 

J C*(k,t)hkC(k,t) d 3 k, 

(201) 

and analogously 

JC*(k,t)^£c(k,t)d 3 k. 

(202) 

Newton's modified second law (|23|) which we have de- 
rived for the non-relativistic case, contains an additional 
"quantum force" Fqp = — VVqp whose expectation 
value equals zero. As a result one arrives at Ehrenfest's 
two theorems. 

< v >= < r >=< V p E{p) > (203) 
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and 



form 



dt <P> 



< F >=< -W > 



(204) 



The salient point here is that these two equations apply 
to the non-relativistic case and we require them to 
persist unaffected in the relativistic case if the particle is 
assumed - as before - to perform a dissipationless motion 
under stochastic extra forces. 

Conversely, one can derive the time-dependent 
Schrodinger equation just by starting from Eqs. (|203|) 
and (|204[ ) and going along the same line of argument used 
in our derivation of the time-dependent Pauli equation 
in Section IXXVI In the following we shall refer to the 
latter. However, instead of 



E(p) = 



2toq 



too c 



V{r) 



we now have 



E(p) = Jp 2 c 2 + to 2 , c 4 +V(r) 



(205) 



= -E fci „+m c 2 



with c denoting the velocity of light in vacuo. 
Hence < V p E(p) >=< v > in Eq. (|203p has to be dealt 
with differently in the relativistic case. Following Dirac 72 
we construct a Fourier-transform H (k) that corresponds 

to the thought-for energy-operator H re i. just as h 2 k 2 in 

-2 

Eq. (|201[) relates to the expression ^ one rewr it es 

Ekin(k) + toq c 2 in Eq. (|205|) in the form 



E k in{k) +m Q c 2 = hc yJE^o fc 2 where p^^hk^ 
and k = 2jp 

and replaces the right-hand side with a 4x4-matrix 
£L Q (k) defined by 

where a denotes constant dimensionless 4x4 matrices, 

the Fourier transform H^ Q (k) must obviously possess the 
property 

3 3 

g^(k) = h 2 c 2 Y, E fc ^ WI = 

p=0 /i'=0 

h 2 c 2 3 3 

— E E *m v % %, + %, «j • 

That means that the matrices have to comply with 
the requirement 

— fa a , + a , a 1 = <L U < 1 . 

As can be verified by just performing the multiplications, 
the matrices meet this requirement if they have the 



and a 




for li = 1, 2, 3 . 



Here a denotes 2 x 2-matrices that are identical with the 

Pauli matrices (|172p . Similar to the latter one can lump 
the 4x4- matrices a together by forming a vector a so 



that H_ Q (k) may be cast as 

H (k) = ca- hk + a Q toq c 2 



(206) 



The feasibility of the above line of thought requires a con- 
sistent extension of the hitherto discussed spinor function 
to a bispinor function 



i>Ur,t) 



where 



^f i) (r,t) = |^,(r,f)|e^m)^, j = l,2. 



1(1) 



The associated phases tp^), Ar,t) represent as in Eq. (|198|) 



potentials of ensemble averages of momenta which means 

<)(^)l 2 



P(r,t)= E 



i=i,2 

(U) 



P(r,t) 



where p^ (r, t) = h V (r, i) . 



,0') 



The quantities pHn are now different for "spin up" 
and "spin down" if the particle in question moves in a 
spatially varying potential. Only in the strictly non- 
relativistic case the spin generating component of the 
quivering motion and the orbital motion remain unaf- 
fected on superposition. In this case we have V'T(I) ( r > = 
|^ m) (r,t)|e^.*). 

If one performs a Fourier transform one obtains in com- 
plete analogy to Eq. (|201|l also in the relativistic case 

<p{t)>= J C_\k 1 t)hkC(k,t)d 3 k. 

Correspondingly one gets 

C?(k,t)Ho(k) C(k,t) d 3 k = 

/ ij)*{r,t) [ca ■ p + a Q m c 2 } ip(r,t) d 3 r . 
■i ~ •• *- ' — 

=Dirac 

We now form < v > according to 
< v >=< V p E(p) >= 

C\k, t) [7T 1 V k Ho{k)} C(k, t) d 3 k . 
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If we substitute Cjk) by its Fourier transform we obtain 
< v >= / ip j '(r,t)caip(r,t)d 3 r . 



Exploiting the identity 

i — r H _,. =— icha. 

^=JJirac ^=JJirac — 

and going through the same set of arguments as with 
deriving the Pauli equation, we arrive at the time- 
dependent Dirac equation 

\Enirac + V ( r M( r ^= ih it( r ^- 



dt ~ 



The derivation can be extended by including electromag- 
netic fields, again in complete analogy to the derivation 
of the Pauli equation. 



XXX. SPATIAL PARTICLE CORRELATION 
BEYOND THE LIMIT OF ENTANGLEMENT. 
SPOOKY ACTION AT A DISTANCE 



As discussed in Section lXXII thc electrons of two hydro- 
gen atoms will respond independently to local perturba- 
tions once the inter-atomic distance has become macro- 
scopically large. The electronic wave function factorizes 
then and becomes the product of two one-particle wave 
functions. One would therefore expect two free fermions 
that have moved sufficiently far away in opposite direc- 
tions with their spins being transverse and anti-parallel, 
to display the same features. If they were still described 
by an anti-symmetric wave function the particle proper- 
ties would remain non-locally intertwined in that each of 
the particles would appear at distant detectors with only 
half of the total probability. Therefore a realistic descrip- 
tion can only be ensured by a product of two one-particle 
wave functions, wavepackets moving in opposite direc- 
tions, one for spin up and the other one for spin down 
or vice versa, the choice randomly distributed among the 
pairs generated in succession. Consequently, there will be 
a complete loss of the "common-cause" -spin correlation 
of the particles when they hit differently oriented spin 
detectors. The latter scatter the incoming fermion de- 
pending on the angle which the fermion's spin direction 
encloses with the scattering plane. To be as concrete as 
possible we refer in this section to the fundamental ex- 
periment by Lamehi-Rachti and MittigS who were able 
to generate pairs of protons of about 8 MeV with spins 
paired anti-parallel and moving apart such that the pro- 
ton's velocities in the center of mass system have the 
same absolute value but opposite directions. The spin 
orientation was analyzed by letting each of the protons 
impinge on a device akin to a Mott detector familiar 
from polarized electron detection. The incoming pro- 
ton is scattered at some carbon atom of a carbon foil. 
Each Mott-type detector is associated with two particle 
detectors whose axes point to the scattering center and 



enclose an angle ±a with the flight direction of the in- 
coming proton. Together with that direction these axes 
form the scattering plane. The differential cross section 
of the carbon scatterer for a proton with spin up perpen- 
dicular to the scattering plane is given by 

a(a,p) = (\f(a)\ 2 + \g(a)\ 2 ) [1 - S(a) sin/3] (207) 

where S(a) represents the Sherman function for car- 
bon/proton scattering, stands for the azimuthal an- 
gle in the plane perpendicular to the proton flight direc- 
tion and f(a) and g(a) denote the scattering and spin- 
flip amplitude. The latter is associated with spin-orbit 
coupling^ which determines also the magnitude of S(a). 
In view of the objective of this article, we wish to em- 
phasize at this point that Eq. (|207|) is a consequence of 
solving the relativistic Pauli equation, and the experi- 
mentally verifiable results that will be discussed below, 
are another objective consequence which is definitely not 
affected by the process of particle detection (the "mea- 
surement" ) . 

In accordance with the notation familiar from EPRB- 
experiments we denote the Mott-type analyzer at the 
end of the left proton track by A and that at the end 
of the right track by B. Furthermore, the detector on 
the right side of the scattering plane will be charac- 
terized by a "+"-sign, that on the left side by a 
sign. The two particle detectors of each Mott-type an- 
alyzer are located at (3 = Ttt/2, and a was set ~ 50 . 
Hence the difference between the respective differential 
cross sections (the "left-right asymmetry" ) is given by 
Act = cr+ - a- = (|/| 2 + \g\ 2 ) 2S. If the spin of the in- 
coming proton encloses an angle A with the normal of 
the scattering plane, the sin- factor in Eq. (|207p becomes 
sin(A =F 7r/2) = =p cos A with (3 = = F 7r /2 denoting the po- 
sitions of the two particle detectors as before. Thus one 
has Act = (|/| 2 + |g| 2 ) 25* cos A. In order to capture the 
general case, we introduce an orthogonal Cartesian co- 
ordinate system whose ir/y-plane is spanned by the two 
proton tracks before they enter the Mott-type analyz- 
ers. The axis of alignment of the proton spins encloses 
in general an angle tp with the z-axis thus introduced. 
The pertinent orientation angles of the scattering planes 
with respect to that z-axis are denoted by 6 and <f> for the 
normals of the A and B-plane, respectively. That means: 
A^4 = 6 — ip and As = 4> ~ f- To make contact to the 
familiar notation, we define a quantity P^iB) through 



'A(B) 



2[|JT 



= p 



MB) 



= S[-£g±\ cosA A(B) ] (208) 



which has the property 



P 



MB) 



p 



MB, 



= 1 



Obviously, P^(b) * s P ro P or tional to the count rate of the 



respective detector, and ^ (Pa(b) 

describes the degree of spin orientation of the incoming 



P A{B)) = C 0SA A{B) 
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proton with respect to the normal of the associated scat- 
tering plane. If the spin of the proton impinging on the 
analyzer at A is parallel to that normal, that is perpendic- 
ular to the associated scattering plane, we have = 
and hence ^ (P4 — P^~) = 1. The joint probability of 
finding the proton pair with one of the protons at A and 
orientation angle = 6 — tp and the other proton at B 
with orientation angle As = 4> — p — 7r is given by 

P joint — (Pa ~ Pa ) C^b — ^b ) = 

P++ + P— - P + - - P-+ (209) 

where P ±± = P± and P ±= F = P| P|. Because of 
the definition (|208[) we have 

P++ + P— + P+- + P-+ = 1 . 

To make sure that the count rates refer definitely to pro- 
ton pairs, the counts associated with P ±±: and P ±=F are 
filtered by coincidence electronics. 

Since for principal reasons one has in general S < (in 
the case under study S w 0.7), P+ - P~ = S cos A 
can never become unity even when the particle enters 
the analyzer with its spin perpendicular to the scatter- 
ing plane, that is when A = 0. It is therefore sug- 
gestive to introduce an 5-independent joint count rate 
P joint — -^iP joint which, on combining Eqs. (|208[) and 
fl2T)9"]) . takes the form 

Pjoint{0, <f>,<p)=- cos(6» - tp) cos(0 - tp) . (210) 

In practice the experiments have been carried out with 
the scattering plane of the B-analyzer lying in the x/y- 
plane, which means = 0. Since all proton pairs are 
prepared 100% polarized, that is with their spins aligned 
parallel and antiparallel with respect to the z-axis, we 
have also tp = so that Eq. (|210p simplifies to 

Pjoint = -cos0, (211) 

and this is in agreement with the experimental results. 
We emphasize again that this equation has been obtained 
by assuming a factorization of the two-proton wave func- 
tion which means that the motion of the "A" -proton 
is controlled only by the potentials specifying the "A"- 
analyzer. There is no influence of the potentials that 
belong to "B". Analogous statements apply to the "B"- 
proton. Hence, for each pair of protons there is no corre- 
lation between their respective "A" and "B"- scattering 
processes. However, it has been the objective of the ex- 
periments, as the authors expressly state, to demonstrate 
that there is such a correlation. Yet in order to prove that 
point, the experiments should have allowed a preparation 
of proton pairs with an axis of spin alignment that en- 
closes an angle tp with the z-axis as originally assumed 
above. According to the established terminology that 
angle has to be regarded as a "hidden variable" . The 
values of tp associated with the various pairs should have 



random character. One can form then a new expression 
from Pjotntid, 4>, <p) by averaging over tp: 

Pav 

(9, 4>) = / p(tp) P Jomt (0, </>, tp) dtp . (212) 

where p(<p) denotes a weight function normalized to 
unity. Clearly, in the experiment the averaging occurs 
automatically and unavoidably. 

If one assumes a uniform distribution of tp over the inter- 
val 7r, that is p(p) — —, and inserts here Pj O i n t(0,<j),tp) 
from Eq. (|210p . one obtains 

P a v(0,<p)=-±cos(6-<P), (213) 

where 

cos(# — tp) cos((f> — p) = \ cos(8 — (f>) + i cos(# + <p — 2<ys) 

has been used. Hence, for = as specified in the 
experiment, Eq. (|213p yields 

Pav(0,(f>) = -| COS0 

which differs from (|21ip by a factor of |. 
At this point it is instructive to contemplate the change 
that would occur if there would be a non-local correlation 
between the two analyzers in the following sense: 
If the "B" -proton has been specified by the "B" -analyzer 
as polarized perpendicular to the associated scattering 
plane, that is if tp> = <fi, and if this property is by some 
"spooky action at a distance" transferred to the "A"- 
proton, tp attains the same value for the "A" -proton. 
The measurement on the "A" -proton would then become 
"contextual" : it would depend on the result obtained for 
the "B" -proton. Consequently, the detection rate (|210p 
would take the form 

Pjoint = — cos(# — 4>) — — cos(a, b) = —a ■ b . (214) 

where we have introduced the quantities a and b as nor- 
mal vectors for the "A"- and "B" -scattering plane, re- 
spectively, which enclose angles 9 and <f> with the z- 
axis. For the situation specified by the experiment (viz. 
4> = 0), this result becomes identical with (|211[) . Thus, 
a distinction between the two mechanisms is not possi- 
ble within the given limitations. One might argue that 
a derivation based on a "spooky- action- at-a-distance"- 
hypothesis has to be rejected anyway. But the same 
hypothesis works perfectly for the analogous experiment 
with pairs of linearly polarized photons where that par- 
ticular limitation does not exist. (S. Aspect et al»2i.) 
By referring to the expectation value 

(^I^a -a®a B - 61*) = -cos(a, b) 

where * denotes the anti-symmetric singlet-state two- 
proton wave function and (Ja/b the spin operators, 
Eq. (|214p is commonly discussed as "the quantum me- 
chanical prediction" for the experiment in question. Con- 
sidering all the details of our analysis it is hard to see 



45 



how this expectation value can have anything to do with 
the experiment except that it happens to yield the same 
— cos(a, b). 

We shortly return to the idea pursued by Lamehi-Rachti 
and Mittig in their paper. In order to exclude the pos- 
sibility that their result might accidentally coincide with 
the prediction of a hidden parameter model, they resort 
to Bell's theorem 75 . It refers to quantities of the type 
P av {a, b) in Eq. (f2T2")) which - according to Eq.flHH - be- 
comes equal to — ^ cos(a, b) if ip is uniformly distributed. 
In general the weight function p((p) will be unknown, and 
hence a complete lack of correlation between the "A"- 
and "B" -scattering processes, as implied by our treat- 
ment, will not show up simply as a numerical correction 
factor of the "correlated result" . Bell could show that in 
performing an EPRB-type experiment one is definitely 
dealing with a non-classical (i. e. non-local) particle cor- 
relation if - irrespective of the form of the weight function 
and irrespective of the kind of hidden variable - the fol- 
lowing inequality is violated: 

\P av (a, b) - P av (a, b')\ < 2 \P av (af, V) + P av {af, b)\ , 

where a, a! ,b, b' denote different analyzer settings. In 
fact, the authors succeeded in verifying this violation, 
but it appears to us, because of the limitations discussed 
above, that this result is not absolutely convincing. 

XXXI. CONCLUDING REMARKS 

In summarizing the essence of quantum mechanics 
Wigner states in a fundamental article 76 under the head- 
line "What is the state vector?": "We recognize ....that 
the state vector is only a shorthand expression of that part 
of information concerning the past of the system which 
is relevant for predicting (as far as possible) the future 
behavior thereof. " 

In our view the most impressive success of quantum me- 
chanics in understanding the stability, composition and 
properties of the building blocks of nature consists in 
predicting the systematic order in the periodic table, 
the phenomenon of chemical valency and the ground- 
state properties of molecules and solids. The state vec- 
tor of these systems, the ground-state wave function 
^( r i, r 2, ■ ■ ■ r/v), is a function of the particle coordi- 
nates r%,r2, ■ ■ -Tn in terms of which their Coulomb in- 
teraction enters the calculation of the system's total en- 
ergy. But for every experimentalist there is no doubt 
that these coordinates are fundamentally inaccessible to 
measurement, and hence cannot possibly be regarded 
as "information gained from measurements" . As is am- 
ply demonstrated by modern ab initio-calculations, the 
wave function allows one to determine the total energy 
as a function of nuclear positions, bonding angles, vi- 
brational frequencies, lattice constants, elastic moduli, 
phonon spectra, saturation magnetizations, electric con- 
ductivities etc.. These quantities are in the spirit of the 



common-sense notion true observables whereas the parti- 
cle coordinates remain definitely hidden parameters. As 
we have repeatedly explained, this applies to the eigenval- 
ues of hermitian operators as well, thus putting a serious 
question mark behind "Kochen-Specker"-typeZ£ and "no- 
go" theorems (s. e. g»^) which are all based on the ex- 
asperatingly artificial assumption that "measurements" 
yield eigenvalues or "probabilities for eigenvalues" . Oc- 
casionally a certain awareness of this puzzling inconsis- 
tency surfaces as in a statement of Wigner'sS: 

"All these are concrete and clearly demonstrated limita- 
tions on the measurability of operators. They should not 
obscure the other, perhaps even more fundamental weak- 
ness of the standard theory, that it postulates the mea- 
surability of operators but does not give directions as to 
how the measurement should be carried out. " 

It is deplorable to notice the impropriety with which 
certain advocates of the Copenhagen school of thought 
dismiss supporters of Nelson's attempt on developing a 
"quantum mechanics without observer" as "stranded en- 
thusiasts" (s. Streater— ), and ironically base their crit- 
icism on the old, actually absurd, arguments how inde- 
terminacy enters the theory through measurement and 
how commuting "observables" correlate with the result 
of simultaneous measurements. All this has been iter- 
ated umpteen times although it is well known to every 
experimentalist that exactly these "measurements" are 
inexecutable altogether. With the same insensitivity to 
reality castigators of the Nelson proponents think it fully 
justified to equate the physics of photon-correlation ex- 
periments with analogous, but actually extremely scarce 
experiments with massive particles. We believe we have 
presented ample evidence that quantum mechanics is in 
detail derivable from classical mechanics plus a modified 
physical vacuum by allowing the latter to undergo en- 
ergy fluctuations. Their action on massive particles is 
calibrated by Planck's constant, and despite their pres- 
ence the conservation of energy (and with free particles: 
the conservation of particle momentum) is ensured on 
average. We hope that the present article can contribute 
to an unbiased reassessment of present-day quantum me- 
chanics concerning these two questions: 1. Which ele- 
ments of the old doctrine are obsolete and dispensable? 
2. Does "measurement" really play a particular role in 
quantum mechanics or is its alleged importance simply a 
misunderstanding? 



XXXII. APPENDIX: DERIVATION OF THE 
NAVIER-STOKES EQUATION 

Given an ensemble of N similarly prepared one- 
particle systems we introduce a transition probability 
P M (r, a, t, At) which denotes the probability of a par- 
ticle being in the elementary volume d 3 o~ around a point 
r + a after a time span At if it has been with certainty at 
point r at time t. This transition probability integrates 



4G 



to unity: 



J P M (r,a,t,At) d?a = l. 



The superscript "M" stands for "Markov process" . 
In terms of this transition probability Einstein's law (|22|) 
on the mean square displacement of a particle under the 
action of stochastic forces may be cast as 



07 <?k P 0, d, t, At) d A a = 5 lk 2v At 



where a — (01, 02, 0-3) . 

Here v = rj/mono denotes the "kinematic viscosity" 
of the embedding system, 77 represents the common 
"dynamical viscosity" and no is the particle density of 
the embedding system. 

The Smoluchowski equation^ describes the temporal 
change of the probability density at r caused by the mo- 
tion of the N independent particles under the influence 
of stochastic forces as a result of which the particles 
perform transitions from previous positions r — a to r 

p(r, t + At) = J p(r - a, t)P M (r - a, a, t, At) d 3 a . 



=G(r-S,S,t,At) 



(215) 



One may evaluate the integral on the right-hand side by 
approximately replacing the integrand G(r — a, a, t, At) 
with a Taylor polynomial of second degree 

G(r - a, a, t, At) = G(r, a, t, At)- 

lEu' 7 ' a kQ^G(r,a,t,At). 
Inserting this expression under the integral of Eq. (|215| ) 
yields p(r, t + At) = p(r, t) - ELi T^W' *) 



x J a k P M (r,a,t,At) d 3 a}- 



-Veh(r,t) At 



iEuaifcWr,*) J a l a k P M (r,a,t,^)d 3 a} 



=6i„ 2uAt 

This may be recast as 
d p 

ttt + V • j c - vAp = where j c = pv c . (216) 
at 

"Fokker-Planck Equation" 

Apart from obeying the Fokker-Planck equation the sys- 
tem of N particles must satisfy the equation of continuity 
as well 



dp 

at 



V • j = where j = pv . 



On forming the difference of these two equations one ob- 
tains 



This is equivalent to 

jd = -^Vp "Fick's Law" . 
If one writes the diffusion current density jd in the form: 

jd — pu (u— "osmotic velocity"), 
one may recast u as in Eq.© 

M(r ' t) = ^^) Vp(r ' t} - 

Thus we have v — v c + u which is just Eq.© used in 
advance in Section HVl 



Replacing p(r,t) in Eq. (|215p with pv c k, one obtains 
dp v ck 



dt 



-\scatter = ~V • (pV c k V c ) + V A(pV ck ) , 



There is an additional (local) change in time of the mo- 
mentum current density effected by the external force 

dP Q t ° k I force = fk(r) = ^-p{r)F k {r) , 
Invoking the Fokker-Planck equation one can write the 



Sum " r Qi" \ force + dP g t ° k \ scatter ™ the form 



dp v c k 
8 

—^ + (v c + 2u)-Vv c -isAv c = — F(r). (217) 
at mo 

Substituting here v c by v — u we arrive at our Eq. : 
^ (v - u) + [(« + u)V(v -u)]-v A(v - u) = i F(r) 

In hydrodynamics \u\ is usually neglected compared to 
\v\ and ^- F(r) is given by the internal mass-referenced 
force ~ p(r,t) Vp(r,t) with p(r, t) — m p(r,t) denot- 
ing the massive density and p(r,t) the pressure. Hence 
Eq.© takes the familiar Navier-Stokes-form: 

p(r, t) + v ■ J -/i Av + Vp (r, t) = 



where p = v p . 
The derivation of this equation is due to Gebelcin 25 . 



V • {jd + v Vp) = where j d = j - j c 
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80 The possibility that both electrons accumulate in one of the 
atoms can safely be excluded. In such a case the other atom 
would be left ionized requiring an energy AE of about 
lRyd. Within a time At that excess energy must disappear 
again where At results from AE At ps h. This yields At ~ 
5 ■ 10~ 17 s, thus excluding the possibility for one of the 
electrons to go back to the ionized atom 10 cm away at a 
speed well below light velocity. 

81 We side here emphatically with John BeU 3 -^ who pleads in 
his article "Against Measurement" for more common sense 
in describing what is actually happening: the time evolu- 
tion of a particular experiment. 

82 "... What we are dissatisfied with is basically not that the 
old perceptions have failed but that they could not be sub- 
stituted by something immediately comprehensible." 

83 However, we want to modify this debatable statement by 
saying that also free electrons display a magnetic moment 
when they are exposed to a magnetic field where their mo- 
tion perpendicular to the field becomes confined to a cir- 
cular area of a certain diameter. 

84 Although spin-orbit coupling represents a constituent of 
the relativistic Pauli approximation to the Dirac equation, 
we consider it here as given. 



